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The  effect  on  rocket  accuracy  of  non-rigidity,  and  especially  of  the 
motor  bending  -when  heated  by  the  propellant  gas,  has  been  investigated 
theoretically  and  experimentally.  The  experiments  were  carried  out  with  a 
Service  rocket  (Rocket  Aircraft  3-inch  No.  1)  which  was  chosen  as  the  most 
suitable  for  the  purposes  of  the  investigation*  It  is  demonstrated  how  this 
rocket  (as  typical)  can  be  modified  so  that  this  source  of  inaccuracy  is 
suppressed.  Charge  designs  in  which  the  propellant  gases  have  free  access 
to  the  motor  body  wall  (Xll),  limited  access  (inhibited  star-section  colloidal 
propellants),  or  no  access  until  late  in  burning  (plastic  propellant)  are 
considered.  Limited  access  may  cause  as  much  distortion  as  free  access. 


The  report  describes  only  mathematical  and  laboratory  work;  the  results 
of  projection  trials  for  dispersion,  which  should  show  the  effect  of  removing 
this  source  of  inaccuracy,  will  be  dealt  with  in  a  subsequent  report. 
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CONFIDENTIAL 

1 .  INTRODUCTION 

Many  attempts  have  been  made  in  the  past  to  calculate  the  ore  t i cally 
the  dispersion  of  unrotated  rockets.  In  these  so  called  tolerance  theories 
it  is  assumed  that  the  rocket  behaves  as  a  rigid  body  and  the  dispersion  is  due 
to  the  measurable  mechanical  tolerances.  The  angular  malalignment  of  the 
nozzle  exit  cone  axis  from  the  rocket  axis  gives  theoretical  trajectory 
deviations  considerably  in  excess  of  those  from  the  other  tolerances.  The 
emergent  gas  jet  derives  almost  the  whole  of  its  momentum  from  the  nozzle  and 
it  can  therefore  be  assumed  that  the  line  of  thrust  is  along  the  nozzle  axis. 

If  this  axis  is  malaligned  so  that  it  does  not  pass  through  the  centre  of 
gravity  of  the  rocket,  there  is  a  turning  moment  on  the  rocket  which  causes 
it  to  turn  from  the  trajectory  i.e.  the  rocket  is  yawed.  The  main  thrust 
then  propels  the  rocket  away  from  the  intended  trajectory.  The  methods  used 
in  the  calculation  of  these  deviations  are  given  in  detail  in  ref.  5*  For 
the  Rocket  Aircraft  3-inch  No.  1 ,  a  nozzle  exit  cone  angular  malalignment  of 
10~3  radians  deviates  the  final  trajectory  0.8  degs.  from  the  line  of  launch. 

The  mean  measured  malalignment  for  a  large  number  of  these  rockets  is 
0.5  x  10**'  radians  approximately  corresponding  to  a  trajectory  deviation  of 
0.4  degs.  This  is  a  circular  distribution  since  the  rockets  are  placed  on 
the  launcher  with  random  orientation  and  it  is  multiplied  by  2/n-  to  reduce  it 
to  a  linear  mean  deviation  (i.e.  deviations  in  one  plane).  Thus  the  expected 
linear  mean  deviation  is  0.3  degs.  However,  the  mean  dispersion  of  this 
rocket  when  filled  tubular  charge  as  obtained  from  a  large  series  of  firings 
carried  out  at  P.D.E.  Aberporth  in  1944  is  1.4  degs.  This  is  appreciably 
greater  than  the  value  0.3  degs.  Similar  discrepancies  are  noted  with  almost 
all  rockets.  There  is  little  reason  to  doubt  the  fundamental  accuracy  of  the 
ballistic  theory  employed.  For  example  the  ratios  of  the  dispersions  of  the 
various  3-inch  rockets  as  calculated  by  this  method  for  the  12  ft.,  24  ft. , 
and  48  ft.  projectors  are  in  reasonably  good  agreement  with  the  observed  values, 
also  the  ratios  of  the  dispersions  with  different  fin  sizes  are  given  correctly. 
The  theory  also  gives  the  correct  wind  constants.  We  can  thus  employ  the 
theory  with  some  confidence  provided  we  know  the  disturbing  forces. 

Three  fundamental  assumptions  are  made  in  the  tolerance  theory. 

1 .  The  rocket  behaves  as  a  rigid  body. 

2.  The  resulting  thrust  vector  arising  from  the  emerging  jet  lies 

along  the  nozzle  exit  cone  axis. 

3.  The  rocket  is  launched  perfectly  so  that  the  angle  of  the  rocket 

axis  and  its  angular  velocity  at  the  moment  of  launch  are 
constant  from  rocket  to  rocket. 

It  is  with  the  first  of  these  assumptions  that  we  are  concerned  in 
this  report.  The  second  and  third  assumptions  are  considered  in  later  reports. 

For  the  practical  investigation  it  is  necessary  to  choose  a  specific 
design  of  rocket  on  which  to  work.  The  design  chosen  is  the  Rocket  Aircraft 
3-inch  No.  1.  This  particular  rocket  has  been  in  Service  since  the  early  days 
of  the  weir  and  it  is  considered  to  be  particularly  useful  for  the  present 
investigation  because  (a)  an  adequate  supply  of  components  is  available,  (b)  it 
has  a  long  history  of  dispersion  and  other  trials  so  that  its  behaviour  is  well 
established,  (c)  a  variety  of  charge  designs  can  be  adapted  to  it,  including 
designs  which  admit  a  copious  flow  of  propellant  gas  against  the  wall  of  the 
motor  body. 

A  detailed  theory  has  been  developed  by  means  of  which  the  motion  of 
a  non-rigid  rocket  is  examined  during  burning.  If  the  motor  tube  bends  during 
flight,  the  line  of  the  nozzle  axis  is  displaced  from  the  centre  of  gravity 
of  the  rocket  thus  producing  a  turning  moment  on  the  rocxet  which  causes  it  to 


CONFIDENTIAL 


1 


yaw  and  deviate  from  the  expected  trajectory.  The  motor  tube  is  under  high 
internal  pressure  and  also  large  set-back  forces  in  flight.  It  can  be  shown 
that  due  to  mechanical  tolerances  mainly  in  the  tube  wall  thickness  these 
forces  bend  the  tube.  Neglecting  temperature  effects  and  taking  the  Young's 
modulus  of  the  steel  to  be  that  at  normal  temperatures,  it  is  shown  that  the 
calculated  bending  due  to  the  measured  mechanical  tolerances  gives  rise  to 
deviations  in  flight  which  are  small  compared  with  deviations  observed  in 
trials.  Thus  the  assumption  of  rigidity  of  the  rocket  tube  is  valid  if  we 
neglect  the  effect  of  temperature.  However,  it  is  found  that  the  asymmetric 
temperature  distributions  around  the  tube  circumference  arising  from  the 
various  factors  discussed  later  give  rise  to  more  serious  bending  of  the  tube 
and  it  is  calculated  that  the  deviations  resulting  from  these  sources  are  of 
the  same  order  as  the  deviations  observed  in  flight.  This  effect  is 
undoubtedly  a  major  source  of  the  dispersion  of  this  rocket. 

There  are  various  factors  which  give  rise  to  the  asymmetric  temperature 
distribution  around  the  tube  section.  To  illustrate  the  orders  of  magnitude-, 
the  following  examples  are  given  for  the  3-inch  rocket. 


X 


(i)  If  the  upper  side  of  the  tube  (i.e.  the  section  above  the  plane  AB) 
is  approximately  12°C  hotter  than  the  lower  side,  the  resultant 
bending  of  the  rocket  is  equivalent  to  offsetting  the  nozzle 
axis  by  10“3  radians  thus  giving  a  mean  circular  deviation  in 
flight  of  0.8°  approximately. 

(ii)  In  the  diagram  above  the  strip  of  breadth  x  ins.  (x  being  small) 

has  a  length  equal  to  that  of  the  tube.  If  the  mean  temperature 
of  this  longitudinal  strip  differs  from  that  of  the  remainder  of 
the  tube  by  an  amount  35/x°C ,  the  resultant  bending  is  again 
equivalent  to  offsetting  the  nozzle  axis  by  10“3  radians  giving 
a  deviation  in  flight  of  0.8°  for  the  3-inch  aircraft  rocket. 

These  temperature  differences  are  small  when  compared  with  the  tempera¬ 
tures  attained  by  the  tube.  Some  typical  values  of  the  latter  for  tubular 
SU/K  charges  for  tubes  with  and  without  a  refractory  insulating  coating  are 
given  in  Table  I,  the  quoted  temperatures  being  those  at  the  middle  section  of 
the  tube. 
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TABLE  I 


Time 

secs. 


Theoretical  temperatures  attained  by  the  walls  of  the 

3-inch  rocket  tube 


Inner  surface 
of  coating 


Temperature  of  coated  tube 

°C. 

Inner  tube  Outer  tube 
surface  surface 


Temperature  of  uncoated  tube 


Inner  tube 
surface 


Outer  tube 
surface 


61 

02 

.  e3 

02 

e3 

0.2 

870 

210 

130 

370 

190 

0.4 

750 

280 

220 

470 

350 

0.6 

710 

340 

280 

550 

44-0 

0.8 

690 

370 

330 

610 

520 

The  effects  are  not  quite  as  serious  as  they  first  appear;  the  bending 
calculated  in  (i)  and  (ii)  is  based  on  the  assumption  that  the  temperature 
differences  are  constant  along  the  whole  length  of  the  tube,  sui  unlikely 
condition  in  practice.  The  principal  factors  which  cause  asymmetric  heating 
of  the  tube  are 

(i)  slight  asymmetries  in  the  refractory  coating  thickness, 

(ii)  charge  asymmetry,  particularly  lateral  displacement,  which  leads 
to  non-uniform  gas  flow  across  the  section, 

(iii)  possible  variations  in  the  roughness  factor  of  the  inside  of  the 
tube, 

(iv)  variations  in  the  wall  thickness  of  the  tube  around  the  section. 

The  tube  coating  consists  of  0.006  ins.  refractory  coating  with  a  top 
protective  bakelite  varnish  coating  of  0.002  ins.  On  the  assumption  that  the 
coating  remains  intact,  it  is  seen  that  the  temperature Qdrop  across  the  coating 
varies  from  approximately  900 °C  at  t  =0.1  secs,  to  300  C  at  0.8  secs.  That 
asymmetric  coating  thickness  affects  the  dispersion  was  shown  by  a  dispersion 
trial.  In  this  trial  the  tipper  half  of  the  tube  was  given  a  coating  of 
thickness  0.01  ins.,  and  the  lower  0.002  ins.;  the  increase  of  dispersion  was 
marked. 

The  thermal  stresses  in  rocket  motor  tubes  have  been  examined  in  some 
considerable  detail  by  Oldroyd  (Ref.  6).  He  has  shown  that  in  the  case  of  the 
Rocket  Aircraft  3-inch  No.  1  filled  tubular  charge,  there  is  appreciable 
plastic  flow  on  the  inside  and  outside  surfaces  of  the  tube  wall  due  to  high 
thermal  stresses.  The  effects  are  maximum  in  the  early  stages  of  burning  when 
the  temperature  gradients  are  maximum.  Some  evidence  for  this  comes  from  the 
bending  tests  for  it  was  noted  with  tubular  charges  that  the  frequency  of 
vibration  as  observed  on  the  bending  records  decreased  markedly  at  0.1  secs, 
after  commencement  of  burning,  and  tended  to  recover  after  0.3  secs.  The 
reduced  vibration  frequency  suggests  reduction  in  the  Young's  Modulus  of  the 
tube  and  this  may  be  expected  to  have  an  effect  on  the  bending.  The  normal 
tubes  are  0.08  ins.  wall  thickness.  Increasing  the  wall  thickness  to  0.18  ins. 
markedly  increases  the  thermal  stresses  due  to  the  increase  in  the  temperature 
gradient.  Comparative  measurements  of  the  bending  of  the  two  sets  of  tubes 
show  that  the  0.08  ins.  tubes  bend  much  more  than  the  thicker  tubes  and  in  fact 
the  ratio  is  that  to  be  expected  if  we  assume  Young’s  Modulus  as  constant. 

The  observed  bending  is  therefore  apparently  due  in  the  main  to  asymmetric 
temperature  distribution.  The  records  of  the  bending  tests  together  with  the 
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measurements  of  the  tubes  before  and  after  firing  are  being  analysed  statisti¬ 
cally  to  determine  the  effect  of  thermal  stresses  and  it  is  hoped  to  give  the 
conclusions  in  a  later  report. 

It  should  be  noted  that  the  bending  observed  in  the  experiments 
described  in  this  report  is  not  directly  applicable  to  flight  conditions.  In 
flight,  the  motor  tube  is  constrained  elastically  only  by  the  shell  and  venturi 
assemblies,  whereas  in  static  firings  it  is  clamped  rigidly  at  the  head  end. 

In  flight  the  natural  frequency  of  vibration  is  calculated  to  be  87  cycles  per 
second  while  in  the  bending  tests  it  is  30  cycles  per  second,  assuming  in  both 
cases  the  value  of  Young's  modulus  to  be  that  at  normal  temperature.  However, 
both  these  values  are  sufficiently  high  to  justify  the  use  of  the  mean  bending. 

In  flight  the  charge  is  subjected  to  an  additional  set-back  force  and, 
in  the  case  of  the  tubular  charge,  this  probably  causes  relative  displacement 
of  the  charge  and  tube  as  the  charge  is  elastically  an  unstable  strut.  This 
displacement  would  have  a  serious  effect  on  the  bending.  Cruciform  charges 
on  the  other  hand  are  forced  against  the  tube  and  this  may  prevent  gas  flow 
between  the  arms  of  the  charge  and  the  tube  so  removing  a  serious  source  of 
bending.  Thus  the  bending  tests  may  not  give  a  trustworthy  comparison  of  the 
two  charges  in  flight. 


2.  EXPERIMENTAL  RESULTS 

2.1  Instrumentation  used  in  the  measurement  of  bending 

The  experimental  arrangement  for  measuring  the  amount  of  bending  of  the 
3-inch  rocket  tube  during  firing  is  shown  in  Figure  30.  The  head  end  of  the 
motor  tube  is  clamped  rigidly  in  two  strong  steel  clamps,  spaced  3-inches  apart, 
the  remainder  of  the  rocket  tube  being  unsupported  so  as  to  permit  free  bending. 
To  determine  the  amount  of  bending  the  lateral  movement  of  the  venturi  end  of 
the  rocket  furthermost  from  the  clamps  is  measured.  For  this,  two  experimental 
techniques  have  been  used,  the  first  based  on  the  use  of  capacity  gauges  and 
the  second  on  magslips.  A  considerable  amount  of  effort  was  expended  in  an 
attempt  to  obtain  reliable  measurements  with  the  capacity  gauges  and  some 
details  are  given  in  order  to  illustrate  the  difficulties  encountered. 

Capacity  gauge  method  Figure  31  shows  the  first  arrangement  used.  A  metal 
ring,  insulated  from  the  rocket  tube  by  a  thin  layer  of  mica  acts  as  one 
electrode.  The  second  electrode  consists  of  a  circular  disc  which  is  fitted 
to  a  coarse  micrometer  screw  thread.  The  condenser  formed  by  these  two 
electrodes  is  placed  in  the  arras  of  a  bridge.  Movement  of  the  rocket  then 
varies  the  capacity  of  the  electrode  system  and  causes  the  balance  of  the  bridge 
to  alter. 


The  electronic  circuitry  associated  with  the  bridge  is  shown  in 
Figure  32.  The  original  unit  was  fed  with  raw  a.c.  to  the  H.T.  line  giving  an 
output  of  rectified  pulses  of  mains  frequency.  Subsequently  the  H.T.  was 
obtained  from  a  stabilised  D.C.  power  supply  giving  an  output  of  radio  frequency 
voltage . 

The  first  valve,  L63  was  an  oscillator  supplying  a  500  Kc.  voltage  to 
the  bridge.  This  was  followed  by  a  two  stage  amplifier  consisting  of  two 
variable  mu  valves  (KTW61).  The  output  from  this  amplifier  was  fed  directly 
to  the  oscilloscope  and  camera.  Movement  of  the  rocket  then  altered  the 
amplitude  of  the  500  Kc.  sine-wave  appearing  on  the  oscilloscope.  Calibration 
was  carried  out  by  screwing  the  electrode  in  steps  of  1/lOth  inch  from  the  ring 
mounted  round  the  rocket.  The  electronic  apparatus  was  situated  in  a  room 
adjacent  to  the  firing  chamber  and  this  particular  type  of  bridge  had  the 
advantage  that  the  screened  leads  to  the  gauge  could  be  made  at  least  30  feet 
long  without  much  decrease  in  sensitivity  or  difficulty  in  obtaining  balance. 
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A  certain  amount  of  success  was  obtained  using  this  apparatus,  in 
particular  the  existence  of  bending  was  established.  The  traces  were, 
however,  susceptible  to  short  violent  kicks  of  a  few  hundredths  of  a  second 
duration  throughout  the  recording.  Although  there  was  no  great  difficulty 
in  reading  the  records  by  extrapolating,  the  presence  of  these  disturbances 
was  unsatisfactory.  The  origin  of  the  disturbances  was  the  rocket  itself 
which  set  up  violent  tremors  in  the  neighbouring  buildings,  shaking  the 
electronic  instruments. 

The  trials  with  the  capacity  gauge  were  exploratory  in  nature  and 
with  the  arrangement  described  above  it  would  not  be  possible  to  measure 
movement  in  two  planes  simultaneously  since  the  gauge  itself  must  be  free  to 
f float'.  Another  capacitance  meter  made  by  Messrs.  Pielden's  Electronics 
employs  a  separate  bridge  unit  and  a  capacity  gauge  with  one  side  earthed. 

It  appeared  that  using  the  arrangement  of  Figure  34,  bending  could  be 
recorded  in  two  axes.  An  initial  difficulty  was  that  beating  occurs  between 
the  two  oscillators  but  this  was  overcome  by  disconnecting  one  of  them  and 
supplying  both  bridge  units  from  one  oscillator.  The  biggest  difficulty, 
however,  was  that  the  lead  from  capacity  gauge  to  bridge  unit  was  only  two 
feet  long  so  that  the  bridges  had  to  be  placed  in  the  actual  firing  chamber. 

The  system  appeared  to  function  in  principle,  but  once  again  the  disturbance 
set  up  by  the  rocket  on  firing  caused  interference.  Extension  of  the  capacity 
gauge  lead  from  the  bridge  unit  decreased  the  sensitivity  correspondingly, 
making  it  necessary  to  raise  the  amplifier  gain  which  in  turn  raised  the 
sensitivity  of  the  recording  apparatus  to  vibration  and  shock. 

Use  of  Magslips 

A  much  more  robust  recording  system  was  sought  and  the  use  of  magslips 
appeared  worthy  of  trial.  These  devices  are  usually  found  as  position 
indicators  in  servo-mechanisms.  It  was  soon  apparent  that  considerable 
advantages  were  associated  with  their  use.  In  particular  the  output- 
displacement  characteristic  was  sensibly  linear  and  their  sensitivity  high, 
so  that  no  auxiliary  amplifier  other  than  that  built  into  the  oscilloscope 
was  needed.  The  operating  frequency  was  1  Kc.  which  enabled  adequate 
frequency  response  (better  than  100  c.p.s.)  to  be  obtained.  Perhaps  the 
biggest  advantage  was  that  calibration  could  be  carried  out  by  movement  of  a 
transmitter  magslip  in  the  recording  room.  The  experimental  arrangement  is 
shown  in  Figure  35«  Typical  records  are  shown  in  Figure  33» 

The  rotor  of  the  transmitter  magslip  is  energised  at  10V  r.m.s.  by  a 
1  Kc.  oscillator.  The  rotor  winding  induces  voltages  into  the  delta  connected 
stator  windings  which  is  passed  on  to  the  receiver  magslip.  Here  the  three 
stator  windings  give  a  resultant  field  so  that  the  voltage  appearing  across 
the  rotor  varies  according  to  its  orientation,  the  output  being  fed  to  the 
oscilloscope. 

The  rocket  is  connected  to  a  lever  system  which  renders  the  magslip 
sensitive  to  horizontal  bending  but  insensitive  to  vertical  bending.  An 
important  feature  is  the  use  of  piano  wire  at  the  right  angle  joint  which 
avoids  the  uncertainty  of  backlash  at  a  pin  joint.  This  useful  device  was 
suggested  by  Mr.  J.  Waddington.  Rotation  of  the  transmitter  rotor  is 
equivalent  to  moving  the  rocket  and  hence  calibration  is  carried  out  by 
movement  of  an  arm  connected  to  the  rotor  to  known  distances,  usually  0.25" 
on  either  side  of  the  zero  line. 

2.2  Discussion  of  the  experimental  results 

In  the  following  pages  are  given  the  results  of  the  bending  experi¬ 
ments  using  both  coated  and  uncoated  3-inch  rocket  tubes.  The  numbers  of 
rounds  in  the  various  categories  fired  are  not  sufficient  to  determine  the 
effect  of  the  coating.  While  it  is  apparent  that  bending  in  the  first  tenth 
of  a  second  is  less  in  the  case  of  coated  tubes,  the  bending  thereafter  in 
each  case  is  too  variable  to  permit  quantitative  comparison.  The  tests  have 
been  carried  out  with  the  various  charges  as  follows 
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Tubular  charges 


Results  of  the  measurement  of  the  bending  of  3-inch  rockets  filled 
with  2.7  -  0.73  inch  tubular  charges  are  shown  in  Figure  4.  The  most  promi¬ 
nent  cause  of  bending  in  this  case  is  probably  the  offsetting  of  the  charge 
from  the  tube  axis.  This  leads  to  several  effects:-  Firstly,  the  rate. of 
burning  of  cordite  in  the  region  of  a  steel  wall  is  reduced  and  hence  this 
wall  is  likely  to  be  the  cooler.  Secondly,  however,  the  conduit  on  the  side 
nearest  the  wall  is  narrower  and  increased  heat  transfer  should  result.  Our 
present  knowledge  of  flow  inside  a  rocket  with  an  offset  charge  does  not 
allow  the  heat  transfer  effect  to  be  calculated,  since  the  flow  is  essentially 
two  dimensional.  We  have  found,  however,  that  rockets  having  charges  offset 
one  tenth  of  an  inch  from  the  tube  axis,  bend  violently  in  the  direction  of 
the  narrow  conduit;  the  results  are  shown  in  Figure  5*  This  offset  is,  of 
course ,  exaggerated  beyond  that  likely  to  be  met  in  practice  though  not  so 
greatly  as  appears  at  first  sight.  There  is  normally  allowed  a  clearance  of 
the  order  0.050  inch  at  70°F. ,  to  ensure  easy  filling  and  reasonable  fitting 
at  the  upper  temperature.  In  flight  the  charge  is  overloaded  as  a  strut  and 
becomes  bowed  between  the  walls  giving  rise  to  conditions  tending  to  those 
obtained  with  the  above  offset  charge  during  static  firing. 

The  tabs  used  to  position  the  charge  are  expected  to  be  the  cause  of 
secondary  effects.  They  are  normally  arranged  in  the  form  of  a  spiral,  each 
tab  causing  a  disturbed  wake  in  the  gas  stream  which  is  not  diametrically 
balanced.  The  rate  of  burning  in  the  region  of  the  tabs  is  altered  by 
erosion  and  variable  stream  pressure.  Solvent  from  the  cement  used  to  attach 
the  tabs  to  the  cordite  causes  a  reduction  in  burning  rate  and  the  cellulose 
acetate  itself  used  for  tabbing  is  attacked  rapidly  by  the  cordite  gases  and 
tends  to  cool  the  stream  locally. 

Opposite  the  tabs  are  small  holes  which  prevent  the  unstable  burning 
which  is  associated  with  secondary  peaks.  In  the  region  slightly  down  stream 
of  these  holes,  temperature  sensitive  paints  show  marked  contours  which  indicate 
that  the  gases  gush  on  to  the  tube  at  these  points.  Figure  6  shows  tempera¬ 
ture  differences  measured  across  a  diameter  10.25  inches  from  the  venturi  end 
of  the  rocket. 

Cruciform  charges 

In  common  with  all  charges  which  expose  the  tube  directly  to  the  main 
gas  flow,  bending  arises  from  wall  thickness  variations  and  insulating  coating 
of  uneven  thickness.  In  addition,  and  characteristic  of  the  cruciform,  is 
the  possibility  of  one  arm  making  poor  contact  with  the  tube,  leading  to  a 
strip  of  the  tube  being  raised  to  a  higher  temperature  than  its  symmetrically 
placed  counterpart.  This  effect  is  difficult  to  avoid  because  of  the  filling 
clearance  needed.  Further,  the  presence  of  different  sized  conduits  must  be 
considered  as  they  may  arise  from  offsetting  of  the  charge  or  as  a  re  suit  of 
extrusion  from  an  uneven  die  during  manufacture.  The  heat  transfer  rate  is 
dependent  on  the  conduit  area.  The  results  of  bending  trials  on  Service 
rockets  using  the  present  standard  XII  charge  are  given  in  Figure  7* 

The  temperature  reached  in  XII  filled  rounds  was  calculated  in  (Ref.l) 
Figure  8.  The  results  of  the  present  thermocouple  measurements  are  given  in 
Figures  9  and  10.  Agreement  with  theory  is  satisfactory,  except  possibly 
near  the  head  end  of  the  rocket  where  the  observed  temperatures  are  somewhat 
higher  than  predicted.  The  positions  of  the  thermocouples  are  shown  in 
Figure  11. 
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Cruciform  charge  in  motor  tube  with  special  insulating  linings 

Firings  were  carried  out  using  tubes  lined  with  0.050  inch  of  an 
insulating  rubber  which  was  fitted  into  the  tube  so  that  the  seam  was  a 
tightly  sprung  joint.  Standard  XII  charges  were  modified  by  planing  0.050 
inch  of  cordite  off  each  arm  before  applying  the  inhibiting  strips  to  allow 
rocan  in  the  motor  tube  for  the  lining.  At  the  end  of  burning  at  least  half 
the  tube  is  still  protected  by  the  lining  the  remainder  of  the  lining  being 
consumed.  The  results  of  this  trial  shown  in  Figure  12  indicate  that  the 
bending  is  much  reduced.  The  rapid  bending  of  one  round  at  the  end  of 
0.6  seconds  is  considered  due  to  the  breakdown  of  insulation  on  one  side. 

An  independent  investigation  of  materials  intended  for  thermal 
insulation  linings  (Ref.  5)  showed  that  polythene  was  very  effective  and  work 
has  been  commenced  to  line  tubes  evenly  with  this  material.  Heating  of  the 
tube  is  avoided  for  about  one  second  by  using  a  0.068"  thickness  of  polythene 
Figures  13  and  14.  A  second  more  readily  available  material  which  showed 
premise  was  rubber  BY16B  (Ref.  5)  s-nd  some  seamless  linings  were  fitted  to 
3-inch  rocket  tubes.  The  material  used  was  poorly  fabricated  and  the  wall 
thickness  varied  from  0.025  -  0.075  inch.  On  firing,  the  tubes  remained 
relatively  straight  until  burn  through  (0.5  sec.)  and  then  bent  violently, 
presumably  due  to  the  insulation  failing  earlier  on  the  thin  side,  Figure  15* 
Temperature  recordings  taken  simultaneously  with  the  bending  trials  confirmed 
that  bending  commenced  at  burn  through. 

Efforts  are  now  being  concentrated  upon  methods  of  producing  an  even 
ooating  of  polythene  for  the  lining. 

Star  centre  charges 

If  no  device  is  used  to  seal  one  or  other  end  of  an  inhibited  charge 
in  a  rocket  then,  as  a  result  of  the  pressure  differential,  gas  flow  occurs 
over  the  inhibiting  coating  on  the  outside  of  the  charge.  This  is  accompanied 
by  heating  of  the  tube ,  characteristically  at  the  head  end  of  the  charge ,  the 
gases  being  cooled  as  they  flow  through  the  narrow  conduit  between  charge  and 
tube  towards  the  venturi.  This  effect  is  most  marked  when  the  conduit  area 
at  the  end  of  the  charge  approaches  that  of  the  venturi  throat  and  the  pressure 
drop  is  high.  Under  these  conditions  the  result  is  occasionally  disastrous  as 
the  inhibiting  coating  is  eroded  away. 

Figure  16  shows  the  distribution  of  maximum  temperatures  in  an  experi¬ 
mental  spin  stabilised  rocket  of  100  mm  calibre  with  open  ended  star  charge. 

The  temperature  rise  at  the  head  end  is  clearly  indicated.  At  the  venturi  end, 
heating  arises  by  conduction  back  fran  the  venturi  plate,  the  temperatures 
indicated  by  the  temperature  sensitive  paints  occurring  several  seconds  after 
the  termination  of  burning.  Similar  results  occur  in  the  3-inch  rocket  when 
filled  with  an  open  ended  star  (Figures  17,  18  and  19).  It  has  been  noticed 
that  the  heating  varies  considerably  in  magnitude  from  round  to  round,  but  the 
temperature  pattern  is  retained. 

In  view  of  the  difficulties  mentioned  above  it  seemed  desirable  to  seal 
the  charge  either  at  the  head  or  venturi  end  so  that  the  differences  in 
temperature  likely  to  arise  around  the  tube  could  be  avoided  or  minimized. 

For  static  firing,  it  is  easiest  to  seal  the  charge  with  a  disc  at  the  head  end; 
although  the  possibility  of  fracture  of  the  charge  on  ignition  is  admitted  in 
so  doing,  no  such  result  has  yet  occurred  in  firings  at  various  temperatures. 

A  trial  using  SU/K  star  centre  charges  resulted  in  very  high  bending 
even  though  the  head  end  of  the  charge  was  sealed  (Figure  20).  Similar  results 
were  obtained  using  am  Improved  charge  in  the  platonised  propellant,  PU 
(Figure  21 ).  These  results  were  surprising  aind  tube  temperatures  were 
measured  in  an  endeavour  to  establish  the  cause.  It  was  at  once  apparent 
(Figure  22)  that  the  temperatures  reached  were  excessive,  and  could  not  be 
explained  purely  on  the  basis  of  gas  flow  into  the  annular  space  around  the 
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charge  followed  by  attainment  of  thermal  equilibrium  with  the  wall  (Appendix  l). 
The  heat  content  of  the  tube  was  some  thirteen  times  too  high  to  account  for 
the  temperature  rise  in  this  way  and  it  appeared  that  excess  gas  flow  must 
have  been  taking  place  towards  the  head  end  during  burning*  The  head— end 
obturator  used  to  seal  the  shell  ring  was  suspected  and  this  was  confirmed  by 
the  fact  that  adhesive  tape  wrapped  round  the  shell  ring  studs  was  set  alight 
on  firing. 

The  3-inch  rocket  was  accordingly  modified  so  that  the  head  end  of  the 
tube  was  sealed  by  a  welded  end  plate.  The  bending  of  the  rocket  was  greatly 
reduced  as  shown  in  Figure  23 •  It  appears  therefore  that  bending  can  be 
largely  prevented  provided  adequate  attention  is  given  to  the  sealing  of  the 
head  end  of  the  rocket  and  the  charge.  The  residual  bending  is  readily 
explained  by  the  fact  that  heating  occurs  from  the  gas  which  equalises  the 
pressure  in  the  conduit  and  the  outer  annulus  between  inhibitor  and  wall. 

Temperature  sensitive  paints  confirmed  that  temperatures  were  low > 
though  conduction  back  from  the  venturi  obscures  the  actual  values  that  arise 
during  burning.  The  temperatures  indicated  by  Termpilac  paints  in  the  cases 
of  ordinary  obturators  and  welded  seals  are  compared  in  Figure  2if. 

Later  trials  with  thermocouples  confirmed  this  view  (Figure  23),  the 
venturi  end  thermocouple  indicated  a  high  temperature  with  this  assembly  as  the 
heat  is  conducted  fran  the  venturi  block.  Service  rockets  must  be 
non-self  propulsive  when  stored  and  the  welded  end  disc  would  not  be  permissible. 
In  consultation  with  E.R.D.E.  the  possibility  of  sealing  the  head  end  obturator 
with  a  cement  was  examined  and  a  polymer  made  from  cashew— nut  oil  and  para¬ 
formaldehyde  CS  which  polymerized  at  air  temperature  was  tried.  Figure  26 
shows  the  results  of  a  bending  using  the  above  mentioned  cement  to  seal  the 
head-end  obturator.  The  bending  is  of  low  order  and  similar  to  that  obtained 
with  welded  head  seals,  but  with  loose  charges  the  bending  is  appreciable. 

Records  are  available  of  star-centre  charges  used  in  the  3-inch  rocket 
during  the  early  part  of  the  late  war  by  the  Low  Pressure  Ballistics  Section 
at  Woolwich. (Ref.  8).  At  the  time,  the  technique  of  inhibition  was  in  an 
early  stage  of  development  and  layers  of  tape,  bonded  with  P.Q.  cement  were 
used.  To  prevent  the  gas  from  ripping  the  inhibiting  tape  away,  a  reversed 
obturator  was  used  at  the  head  end  and  this  appears  to  have  functioned  as  a 
head  end  seal.  Projections  of  these  charges  showed  that  they  had  a  relatively 

low  dispersion;  particularly  so,  because  the  thrust  was  lower  than  that  of 
tubular  filled  rockets  and  of  longer  duration,  giving  a  lower  launch  velocity 
and  raising  the  expected  dispersion.  Sane  of  the  improvement  was  ascribed  to 
the  venturis  which  were  machined  from  the  solid  and  which  gave  a  slightly 
improved  dispersion  with  tubular  charges.  Nevertheless,  several  trials  were 
carried  out  showing  that  stair  centre  charges  were  associated  with  a  low 
dispersion  in  this  particular  rocket. 

Preliminary  work  with  rockets  filled  with  plastic  propellant  is 
described  in  (Ref.  9).  In  this  case  emphasis  is  placed  on  the  centralizing 
•f  the  central  conduit  in  order  to  avoid  extremely  large  end  deflection  of  the 
order  1-2  inches  which  occur  due  to  uneven  sliver  left  at  burn  out.  Lift 
forces  eire  then  liable  to  deviate  the  rocket  since  the  velocity  is  high.  No 
bending  occurs  with  case  bonded  plastic  propellant  filled  rockets  until  the 
sliver  period  since  the  tube  is  completely  insulated. 

3.  DERIVATION  OF  THE  EQUATIONS  OF  MOTION  OF  THE  NON-RIGID  ROCKET  IN  FLIGHT 

The  motion  of  the  rigid  rocket  in  flight  has  been  described  in  detail 
by  many  writers.  The  deviations  of  the  rocket  from  the  flight  trajectory  due 
to  such  factors  as  tolerances  and  wind  have  been  calculated  for  many  weapons. 

An  account  of  the  work  and  the  methods  of  calculation  are  given  by  Davies 
(Ref.  7).  In  the  calculations  given  there  it  has  been  assumed  that  the  rocket 
behaves  as  a  rigid  body  i.e.  it  does  not  bend.  In  this  section  we  extend  the 
theory  to  the  non-rigid  rocket,  i.e.  a  rocket  which  bends  uuder  the  stresses 
arising  and  the  temperature  effects.  The  derivation  of  the  equations  of  motion 


8 


is  described  in  detail  so  that  the  reader  may  comprehend  the  extent  of  the 
analysis  and  the  assumptions  made*  The  method  of  solution  of  these  equations 
is  complex  and  exceedingly  lengthy  and  is  not  given  here;  the  conclusions 
from  the  analysis  only  are  given.  A  detailed  account  of  the  method  of 
solution  will  be  given  in  a  later  report. 


Fig.  3.1 


The  x  axis  is  taken  along  the  central  axis  of  the  rocket  and  x  =  o  is 
taken  at  the  line  of  intersection  of  the  tube  and  shell  assembly.  The  line 
of  intersection  of  the  nozzle  assembly  and  the  tube  we  take  as  x  =  l.  The  y 
axis  is  taken  at  right  angles  to  the  x  axis  and  the  (x  ,y )  plane  is  the  plane 
of  motion.  We  adopt  the  convention  that  the  shearing  force  S  at  any  section 
*  is  the  sum  of  the  external  forces  in  the  vertical  direction  y  on  the  right 
v>nr>H  side  of  x  while  the  bending  moment  B  is  the  sum  of  the  anti-clockwise 
moments  of  these  forces  about  x.  The  wall  thickness  of  the  tube  is  t  and  the 
inside  and  outside  radii  are  r1  and  r2  respectively.  The  shell  assembly  of 
mass  W1  and  the  nozzle  assembly  of  mass  W2  we  treat  as  rigid  bodies. 


The  coordinates  (x  ,y)  specify  the  position  of  an  elementary  length  of 
the  rocket  axis.  In  the  first  instance  let  us  suppose  that  there  are  no 
external  forces  on  the  system  (except,  of  course,  d'Alembert  forces).  We 
suppose  the  rocket  executes  small  motions  about  its  initial  direction.  The 


terns  are  all  taken  as  small  and  we  neglect  their  squares  and  products. 

0X  0  x 


Inside  the  tube  is  the  rocket  charge  which  rests  on  the  grid  at  the  nozzle  end. 
It  is  not  anchored  at  the  head  end  but  plastic  tabs  placed  in  contact  with  the 
tube  and  charge  restrict  their  independent  motion.  In  this  report  we  suppose 
the  tube  and  charge  bend  and  vibrate  in  phase  with  equal  amplitude.  The 
effect  of  this  assumption  is  considered  in  a  later  sypplementory  report. 


The  stress  systems  arising  in  the  rocket  tube  are 


(i)  the  elastic  forces  due  to  the  bending  of  the  rocket, 

(ii)  the  reactions  of  the  gas  flowing  in  the  conduit, 

(iii)  the  inertia  forces, 

(iv)  thermal  stresses. 
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Pig.  3a 


Pig.  3b 


The  section  of  the  tube  in  Figure  3a  is  perpendicular  to  the  x  axis, 
anri  AB  is  an  arbitrary  plane.  The  perpendicular  distances  from  this  plane 
are  denoted  by  z.  The  elementary  length  of  the  rocket  tube  of  which 
Figure  3a  is  the  section  is  given  in  Figure  3b.  In  its  initial  state  this 
section  has  a  radius  Rq  and  subtends  an  angle  0O  at  the  centre  of  the  circle 
of  radius  R.  In  its  strained  state  let  these  values  be  R  and  0  respectively. 
We  suppose  the  temperatures  of  the  elements  of  the  tube  are  functions  of  z  so 
that  we  consider  the  two  dimensional  problem  only.  The  temperature  of  the 
tube  in  the  -initial  state  is  taken  as  the  zero  of  temperature.  If  kz  is  the 

mean  coefficient  of  thermal  expansion  between  the  zero  and  Tz,  the  extension 

of  the  longitudinal  element  at  z,  which  gives  rise  to  stress,  is 

jjCR  +  z)e  -  (kzTz  +  l)(Ro  +  z)e0J/(Ro  +  z)0o  ,  (3.1) 

Tz  being  the  temperature  at  z. 

If  at  a  plane  z  =  a,  the  extension  is  zero 

(R  +  a)0  =  (Ho  +  a)0o(l  +  k&  T&) 

In  terms  of  this,  the  extension  can  be  written 

U  -  a)(i  -  ±  )  .  (ka  Ta  -  kz  Ta)  , 


(3-2) 

(3-3) 


neglecting  the  squares  and  products  of  l/R  and  1/Rq .  If  Ez  is  longitudinal 
Young's  modulus  at  z,  the  force  on  the  element  is  thus 


E 


(z  -  a)(i-i  )  + 


k  T 
a  a 


-  k 


T*j 


dA 


(3.4) 
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dLA  being  the  cross  sectional  area  of  the  element.  The  total  force  across  the 
section  is 

Fx  =  j  \ 

and  the  moment  of  the  forces  about  z  =  z,  is 

j  Ez  [(*  -  a4  -  V  +  ka  Ta  '  ks  Tz] 


[(•- ‘>(4 -$,)+*W  (5<5) 


which  may  he  written 

-  F*z,  -  ^  )  +  T* ,  (3«6) 

where 

li  =  ^  (z  -  a)  EzzdA  ,  T'=  j  Ezz  (k&  T&  -  kz  Tz)  <3A,  (3*7) 

the  integrals  being  taken  around  the  complete  cross-section  of  the  tube* 

The  moment  of  the  forces  about  z  =  z1  ,  is  equivalent  to  a  force  Fx 
parallel  to  the  axis  at  z  =  o  and  a  couple  of  magnitude 

(3-8) 

The  reaction  on  the  tube  of  the  gas  in  the  conduit  can  he  determined 
in  the  following  manner.  Let  us  assume  the  area  of  the  gas  conduit  along 
the  length  of  the  motor  is  constant;  departures  from  this  assumption  can  be 
shown  to  be  of  little  importance . 


Let  p  be  the  pressure  at  CD  in  Figure  3a,  the  pressures  at  GF  is 

p  +  9£  di .  In  an  interval  6t,  the  momentum  crossing  CD  is  pv2  A6t,  p  being 
3i 

the  density  and  v  the  velocity  of  the  gas  stream  and  A  the  area  of  the  gas 
conduit.  In  the  interval  6t  the  momentum  flowing  across  FG  is 

A(  pv2  +  6x  )  6t 

OX 

Now  suppose  the  section  of  the  motor  is  at  rest  and  also  the  gas  is 
liberated  from  the  charge  with  no  momentum.  If  AS  is  the  shearing  force 

exerted  on  the  motor  by  the  gas  stream,  from  the  conservation  of  momentum  in 
the  vertical  plane  we  have 


[pff*-pff*-£<*!f>*-*sP> 

■  ['"*  If  *  &T  "  <"*  If  j  m  ‘ 


or 
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From  the  conservation  of  momentum  in  the  horizontal  direction  we  have 


p  +  pv2  =  P0  , 

P0  being  the  head  end  pressure,  and  thus 


In  this  we  have  assumed  the  section  of  the  motor  to  be  at  rest.  If 
we  include  the  motion  in  the  vertical  plane  we  have  to  add  to  the  above  the 
shearing  force  due  to  the  inertia  acceleration  in  the  vertical  plane  and  also 
the  component  due  to  the  change  in  the  direction  of  the  gas  stream.  It  is 
easily  seen  that 


(3-9) 


In  these  we  have  neglected  the  effect  of  the  forward  velocity  of  the 
rocket.  This  will  be  discussed  in  a  later  report. 

If  w  is  the  weight  of  the  motor  tube  per  unit  length,  the  longitudinal 
force  at  the  section  s  is 


Fx  =  P0 A  -  Wjf  -  wxf 

In  deriving  the  elastic  stresses  we  used  an  arbitrary  reference  plane 
for  z  =  o.  We  have  also  used  a  reference  axis  for  y.  Let  the  z  =  o  plane 
be  specified  by  Y  when  referred  to  the  y  =  o  axis.  The  ordinate  y  is  on  the 
resultant  momentum  vector  of  the  gas  flow  in  the  conduit. 

Consider  the  forces  and  couples  in  the  elementary  length  of  the 

rocket  given  in  Figure  3t>«  On  the  section  CD  we  have  a  vertical  force  -S 

a  clockwise  couple  B,  a  horizontal  force  -Fx  at  z  =  o.  On  FG-  we  have  a 

3S  3B 

vertical  force  S  +  —  dx ,  an  anti-clockwise  couple  B  +  —  dx  at  z  =  o.  On 

3i  3x 

the  element  there  is  an  equivalent  vertical  force  of 


3S„ 

P  dx 


3x 


_3_ 

3t 


dx  • 


Equating  vertical  forces  we  find 


(3.10) 


Taking  moments  we  find 


c  +  9B  _  F  91 

S  +  3*  Px  3x 


=  Ii 


33y 

3x3t2 


+  wfe 


(3.11) 


where  I3dx  is  the  moment  of  inertia  of  the  element  about  an  axis  through  its 
centre  of  gravity  and  perpendicular  to  the  rocket  axis  and  B  is  given  by 
(equation  3*8). 


B  =  -  p(^  -  ^ )  -  T', 


(3.12) 


12, 


with 


1  _  9fZ 

r  “  ax2 


The  term  e  is  the  distance  of  the  mechanical  centroid  of  the  section  above 

z  =  o.  Differentiating  equation  3»1'l  with  respect  to  x  we  obtain  ^  which 

can  be  substituted  into  equation  3, 10  to  obtain  the  differential  equation  of 
motion. 

In  the  preceding  analysis  two  planes  of  reference  have  been  employed 
namely  y  =  o  and  z  =  o,  the  ordinates  of  the  latter  when  referred  to  y  =  o 
being  denoted  by  Y.  Only  the  differentials  of  y  and  Y  occur.  The  y  plane 
is  such  that  at  the  point  x  the  resultant  momentum  vector  of  the  gas  stream 
in  the  conduit  passes  through  y.  The  choice  of  Y  is  still  arbitrary  and 
this  is  chosen  from  later  considerations.  We  have  assumed  in  the  analysis 
that  the  momentum  vector  of  the  gas  stream  is  determined  entirely  by  the 
elementary  section  of  the  motor  so  that  y  expressed  the  position  of  the 
elementary  section  of  the  motor  tube. 

We  have  now  to  derive  the  boundary  conditions.  Let  us  first  consider 
the  contribution  to  the  shearing  force  and  the  bending  moment  from  the  shell 
and  nozzle  assemblies.  W1  is  the  mass  of  the  shell  assembly  and  I1  its 
moment  of  inertia  about  an  axis  through  the  centre  of  gravity  of  the  shell 
assembly  and  perpendicular  to  the  central  axis.  W2  is  the  mass  of  the  nozzle 
and  I2  its  moment  of  inertia  about  an  axis  through  the  centre  of  gravity  of 
the  nozzle  assembly  and  perpendicular  to  the  central  axis.  Let  l 5  be  the 

distance  of  the  centre  of  gravity  of  the  shell  assembly  from  x  =  o  and  l2  the 

distance  of  the  centre  of  gravity  of  the  nozzle  assembly  from  x  =  l.  The 

velocity  of  the  centre  of  gravity  of  the  nozzle  is 


Thus  the  contribution  to  the  shearing  force  at  x  =  l  of  the  inertia  of  the 
nozzle  is 

s2  =  -  w2  jjy  +  i2  ,  *  =  1  (3.13) 


while  the  contribution  to  the  bending  moment  is 


B2  =  W2  I2 


+ 


(3.14) 


Similarly  we  find 
and  bending  moment  of  the 


at  x  =0,  the  contributions  to  the  shearing  force 
shell  inertia  are 


(3.15) 


Let  the  resultant  force  on  the  nozzle 
charge  be  Q  in  the  direction  of  x  increasing, 
the  charge  and  tube  at  x  =  l  is 


sZjl* 

9x01? 


*  =  o. 


(3.16) 


y  the  grid  and  the  base  of  the 
The  total  longitudinal  force  on 
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F  =  P0  A  -  W,f  -  wlf  =  Q  +  W2f . 

•  •  i 

This  is  the  equation  for  f  the  acceleration.  Owing  to  the  curvature 

of  the  tube  we  have  an  additional  shearing  force  at  s  =  t  of  (Q  -|^)  ; 

91  x  =  l 

this  we  add  to  S2.  Similarly  we  add  a  shearing  force  of  P0A(|^)  to  S1 . 

1=0 

Let  us  suppose  the  centre  of  pressure  at  the  head  end  is  at  a  distance 
et  above  z  =  o,  thus  we  add  to  Bi 

-  w,f  l,  (§)  -  P0Ae2  +  W^Ce,  +  e2 ) 

t=o 

where  e3  is  the  distance  of  the  C.  of  G.  of  the  shell  assembly  above  the  centre 
of  pressure  in  the  static  undeformed  state.  This  is  made  more  definite  when 
the  axes  have  been  defined.  Similarly  to  B2  we  add 

-  W2fl2  (§J)  -  Qe4  -  W2f  (e5  +  e4)  , 

I  =  l 

where  e4  is  the  distance  of  Q  above  the  plane  z  =  o  and  e5  the  distance  of  the 

centre  of  gravity  of  the  nozzle  assembly  above  Q.  The  exact  specification  of 
the  terms  is  given  later  when  the  axes  Eire  defined. 

Further  if  due  to  the  asymmetry  of  the  nozzle ,  the  line  of  thrust  is 
inclined  at  an  angle  +  and  displaced  laterally  upwards  above  z  =  o  by  an 
amount  e4  we  add  -  F^  to  S2  and  Fe6  to  B2 

Finally  the  gas  flow  in  the  nozzle  gives  a  shearing  force  at  i  =  l  of 

-  2Apvl3 


and  a  bending  moment  of 

“  lift  * 

l3  being  the  length  of  the  nozzle.  These  expressions  we  add  to  S2  and  B2 
respectively. 

In  addition  to  these  we  have  external  forces  and  couple.  We  denote 
these  by  a  lateral  body  force  of  j(x  ,t)  per  unit  length  and  latereil  forces 
and  couples  of  K^t),  H^t)  amd  K^t),  HaCt)  at  x  =  o  and  i  =  l  respectively. 

The  resulting  values  of  S1 ,  B1  and  S2 ,  B2  we  put  equal  to  the  values 
of  S  and  B  given  by  equations  3*13  onwards  to  obtain  our  boundary  conditions. 

To  illustrate  the  form  of  the  solution  the  differential  equation  and 
boundary  conditions  are  derived  for  a  perfect  rocket.  The  complete  specifi¬ 
cation  and  solution  of  the  differentisil  equation  and  boundary  conditions  are 
given  in  a  later  report. 

The  following  assumptions  were  made. 

(i)  The  mechanical  centroid  of  any  section  coincides  with  the  centre 
of  pressure  in  the  conduit, 

(ii)  The  Young's  modulus  is  constant  around  the  section. 
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In  this  case  we  take  the  reference  axix  z  =  o  as  the  meohanical 
centroid  of  the  section  so  that  the  axis  of  y  and  Y  coincide*  Also  we  have 


/ 


z  Ez  dA  =  0, 


H-zj 


z2  dA 


and 


T  =  -  E  f  kz  Ts 


z  dA  , 


thus  if  we  assume  the  conductivity  and  the  temperature  to  be  completely 
symmetrical  about  the  centroid  of  the  section  we  have  T  =  0.  Also  we  omit  ’v 
the  "gas  whip”  terms  in  Apv. 

With  these  assumptions  the  differential  equations  and  boundary 
conditions  become 

nD*y  +  (w*  +  wt)fD2y  +  wfDy  -  l3  D2y  +  w  =  (*»*)»  o  <  x  <  l  (3«17) 

-  /iD3y  -  WifDy  +  I3  Dy  =  Wi  jjj^r  “  li  ^2  DyJ  “  *  =  o,  (3«1 8) 

-  n&y  +  W2fDy  +  I3  Dy  =  ~  W2  f|^  +  l2  Dyj  +  ^(t),  x  =  l,  (3-19) 


v&y  =  -  w2 12 


\iT?y  =  «  W1  li 


Dy 


Dy 


+  ii 


32 

at2 


Dy 


W2fl2  Dy  +  He(t),  x  =  l, 
By  ~  HfCt),  x  =  0. 


(3.20) 


(3.21) 


The  equations  of  motion  have  been  solved  for  the  Rocket  Aircraft 
3-inch  No*  1  rocket  with  parallel  sided  shell  and  A*D*  fuze  fired  from  a  two 
railed  launcher  12  ft.  in  length.  In  the  numerical  calculations  the  value  of 
Young's  modulus  has  been  taken  to  be  that  for  a  cold  tube  so  that  the  heating 
of  the  tube  has  been  neglected.  The  conclusions  from  the  analysis  are  as 
follows 

(i)  The  fundamental  frequency  of  the  motor  tube  is  87  cycles  per 
second. 


(ii)  As  a  result  of  this  high  fundamental  frequency  the  deviation  of 
the  rocket  due  to  small  lateral  forces  is  almost  identical  to 
that  of  a  rigid  rocket. 


(iii)  The  bending  due  to  the  high  internal  stresses  and  the  tolerances 
produce  deviations  which  are  small  compared  with  the  deviations 
in  flight;  the  maximum  mean  deviation  calculated  being 
0.03  degs. ,  this  being  due  to  the  known  wall  thickness 
variations. 


(iv)  The  velocity  at  which  the  tube  "flutters"  due  to  the  aerodynamic 
forces  is  considerably  in  excess  of  the  velocities  attained  by 
this  rocket.  (At  the  flutter  velocity  resonance  vibration 
takes  place). 


15 


In  the  calculations  it  has  been  assumed  that  the  fins  are  rigid,  no 
account  has  been  taken  of  fin  vibration. 

Finally  we  conclude  from  our  analysis  that  if  the  heat  is  kept  from  the 
tube,  the  assumption  of  non— rigidity  made  in  the  classical  rocket  theory  is 
valid.  In  the  next  section,  however,  it  is  shown  that  the  temperature  effects 
are  serious. 

Section  4.  The  bending  of  the  rocket  due  to  temperature  effects 

In  the  previous  section  we  have,  in  the  numerical  analysis,  assumed 
that  the  temperature  is  constant  around  a  section  of  the  motor  and  also  that 
Young's  modulus  is  constant.  In  view  of  the  high  temperatures  attained  by  the 
tube  we  now  investigate  temperature  effects. 


We  have  seen  that  with  a  perfect  rocket  i.e.  uniform  wall  thickness 
and  symmetrical  temperature  distribution  the  pressure  has  no  effect  on  the 
bending.  Further  the  set-back  forces  create  an  effective  compression  on  the 
tube  but  with  the  3-inch  rocket  already  investigated  they  have  little  effect  on 
the  motion. 


In  the  first  instance  we  assume  the  rocket  to  have  a  uniform  wall 
thickness  but  that  the  tube  has  a  temperature  variation  around  the  section, 
there  being  no  variation  along  the  length.  We  suppose  the  temperature 
variation  stiff iciently  small  to  have  a  negligible  effect  on  the  Young's  modulus. 


As  before  we  take  the  z  =  o  axis  such  that 


dA  =  0  so  that 


dA  , 


T*=  -  k  E 


TzdA  . 


The  bending  moment  at  the  section  is 


B  =  (i/R  +  t'  , 

and  if  we  neglect  the  effect  of  the  set  back  forces  the  mean  bending  couple  B 
is  zero  so  that 

1/R  =  -  t'/u  . 


\  *  / 
\  / 


Fig.  4.1 


Fig.  4.2 
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Assuming  the  tube  to  be  circular,  r  the  mean  radius  and  a  the  mil 
thickness  we  have 

r 

1/R  =  k  I  T  r2  sin  0  a  d  0  /  zir 

Jn 


Jar 3  sin2  0  d  0  , 


2ir 

1 


T  sin  0  d  0  . 


Let  us  suppose  the  upper  half  of  the  tube  is  T  degrees  above  the  lower 
half  i.e.  the  temperature  above  A  B  is  T  degrees  higher  than  that  below. 

Hence 


v 


1/R  =  ^  J  sin  0  d  0  =  2  k  T/n-  r  . 
»o 

The  angle  0  subtended  by  the  rocket  motor  tube  is 


0  =  l/R  =  2'k  l  T/irr  . 


For  the  3-inch  rocket  with  parallel  sided  shell  and  A  D  fuze  it  can  be 
shown  that  if  the  angle  subtended  by  the  tube  is  0,  the  centre  of  gravity  is  at 
a  perpendicular  distance  from  the  line  of  thrust  which  to  a  very  good  approxi¬ 
mation  is  l  0/2.  Thus  the  effective  nozzle  offset  is  0/2,  where 

0/2  =  k  l  T/ir  r  (4.1) 

Taking  k  =  10**^,  l  =  3*6  ft.,  r  =  1.59  ins.  we  find 

0/2  =  8.64  T  x  icf5. 

An  effective  nozzle  offset  of  10*"^  radians  causes  a  deviation  in  flight 
of  0.8  degs.  Thus  the  deviation  in  flight  is  a  where 

a  =  8.64  x  0.8  x  T  x  10  2, 

and  thus  a  temperature  difference  of  14.5°C  between  the  two  halves  of  the  tube 
causes  a  deviation  in  flight  of  one  degree. 

If  a  strip  of  width  i  along  the  whole  length  of  the  rocket  tube  were 
at  a  temperature  T  above  or  below  that  of  the  remainder  of  the  tube,  the 
subtended  angle  is 

0  =  l  k  T  x/ir  r2 


assuming  x  is  sufficiently  small  for  sin  0  to  be  taken  as  unity.  The  effec¬ 
tive  nozzle  offset  is 


0/2  =  k  l  T  x/2  7r  r2 
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and  the  deviation  in  flight  is 

a  =  0.8  x  10"2  l  T  x/2  ir  r2  (4.2) 

Taking  i  =  1  inch,  a  temperature  difference  of  46°C  causes  a  deviation 
in  flight  of  one  degree. 

These  temperatures  appear  small  when  compared  with  the  high  tempera¬ 
tures  attained  in  the  tubes  of  the  3-inch  motors  but  it  must  be  noted  that  we 
have  stipulated  these  temperature  differences  as  applying  along  the  whole 
length  of  the  tube,  an  unlikely  condition  in  practice. 

Measurements  of  the  wall  thicknesses  of  a  large  number  of  tubes  have 
shown  that  the  wall  thickness  varies  systematically  around  the  section  and  this 
variation  is  maintained  with  only  slight  change  along  the  whole  length  of  the 
tube.  The  results  of  the  measurements  are  given  in  appendix  I.  If  with 
tubes  of  this  type  we  chose  the  reference  plane  z  =  o  to  coincide  with  the  line 
of  centroids,  the  centre  of  pressure  in  the  conduit  no  longer  coincides  with 
z  =  o.  The  rocket  is  bent  by  two  factors 

(i)  the  temperature  variation  resulting  from  the  uneven  wall  thickness, 

(ii)  due  to  the  separation  of  the  centre  of  pressure  from  z  =  o,  the 
gas  pressure  exerts  a  bending  couple  on  the  tube. 

With  uncoated  tubes  and  tubes  coated  with  0.006  ins.  of  refractory 
coating,  the  former  is  markedly  greater  than  the  latter  and  for  the  moment  we 
concern  ourselves  with  the  effect  of  the  temperature  variation. 

We  choose  our  reference  plane  for  z  such  that 

j  z  Ez  dA  =  0  ,  (4.3) 

and  this  defines  our  reference  plane  uniquely.  Thus  we  find 


U 


Ez  dA  , 


z  k  T  dA  . 
z  z 


(4.4) 


Let  us  suppose  the  heat  transmission  coefficient  from  the  gas  to  the 
metal  is  constant  all  around  the  section.  During  the  time  of  burning  there  is 
negligible  heat  flow  around  the  tube,  the  heat  flow  being  almost  entirely  radial 
and  the  local  mean  temperature  of  the  tube  is  dependent  on  the  wall  thickness 
at  the  point.  For  example  in  Figure  2a,  the  mean  temperature  of  the  element 
dA  is  dependent  on  the  wall  thickness  at  that  point. 

The  temperature  depends  on  the  depth  u  from  the  inside  surface.  The 
value  of  z  is  practically  constant  over  the  element  and  by  the  temperature  Tz 
we  imply  the  mean  temperature  over  dA  or 


(4.5) 


L  being  the  local  wall  thickness. 

Let  L  be  the  mean  wall  thickness  over  the  cross-section  of  the  tube  and 
let  0m  (l)  be  the  mean  wall  temperature  with  this  value  L.  Thus  for  small 

variations  of  L  about  L 
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r 


T  d  us  T  L 
z 


=  L  [em  (L)  ♦  15*  (L  -  L)J  . 


(U6) 


It  is  shown  from  the  measurements  of  wall  thicknesses  of  a  large 
number  of  tubes  given  in  the  appendix,  that  the  wall  thickness  at  z  is  given 
by  a  relation  of  the  form 


L  =  L  +  b  (z  -  z2)/r 


(4.7) 


r  being  the  radius  of  the  tube. 
Thus  we  have 


z  Tz  dA  =  6m  (L)  j)  z  d  A  +  £^m  |j-  j)  (z  -  z2  )  zdA 


■0 


6m  (L)  -  — ^2-  ^^m^  |zdA  +  5.  j)  z2dA 


(4.8) 


The  Young*  s  modulus  Ez  is  approximately  a  linear  function  of  tempera¬ 
ture  and  we  write 


Ez  =  Eq  +  E1  Tz  , 


(4.9) 


and  thus 


f 


EzzdA=E^^zdA  +  E1  |  z  Tz  d  A 


•r 


j  Ez  z  d  A  =  Je,,  +  E1  j^6m  (L)  -  |  z2  fd6' 


&}}} 


z  d  A 


+  E. 


90  b  I 

f 


z2  d  A  , 


(4.10) 


and  thus  since  the  left  hand  side  has  been  equated  to  zero  to  specify  the  axis 
we  have 


f 


z2  d  A 


z  d  A  =  —  9L 


•m  i  <S 

.  r  J _ 

Eo  +  (L)-  —  z2 


(4.11) 


The  values  of  E/Eq  ,  the  ratio  of  Young's  modulus  at  the  specified 
temperature  to  that  at  zero,  together  with  the  values  of  k  and  k  E/E 0  for 
medium  carbon  steel  are  given  in  the  following  table. 
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?emp.  (°C.) 

E/Eo 

k  x  105 

kE/E„  x 

0 

1.0 

1.15 

1.15 

50 

0.985 

1.17 

1.18 

100 

0.970 

1.20 

1.21 

150 

0.955 

1.23 

1.24 

200 

0.94 

1.25 

1.27 

250 

O.92 

1.27 

1.29 

300 

0.90 

1.30 

1.30 

350 

0.88 

1.33 

1.32 

400 

0.86 

1.33 

1.33 

450 

O.83 

1.37 

1.33 

500 

0.81 

1.40 

1.32 

550 

0.71 

1.43 

1.21 

600 

_  0.60 

1.45 

1.05 

The  product  kE  yaries  with  temperature  and  the  expression  for  T  given 
by  equation  4*4  becomes  complicated.  It  is  considerably  simplified  if  we 
suppose  kzEz  is  a  function  Tz  as  given  by  the  above  table.  Let  us  write 


E  k  =  (Ek)  +  I—  (T  -  6  )  , 
z  z  '  'm  3T  '  z  nr  * 


(4.12) 


(Ek)_  being  the  value  at  0  .  The  moment  of  area  I  of  the  section  about  z  =  o 
•  in  m 

is 


-f 


z2dA  , 


and  in  equation  4.11  the  term  bz2 -  can  be  neglected  compared  with  0  and  thus 

9L  m 


f 


z  d  A  =  -  E<  bl  /E  r 
9L~*  7  m 


(4.13) 


Em  being  the  Young's  modulus  at  the  temperature  ©m  (L). 


In  appendix  2,  it  is  shown  that  for  the  3-inch  rocket  tube,  wall 
thickness  0.08  ins. 


90 

C-J£) 

9L  L=L 


=  -  11.2  (0,  +  0a)/2  , 


01  and  02  being  the  inner  at  outer  temperature  of  the  tube  wall.  Thus 


f 


z  d  A  =  5*6  (01  +  02)E1  b  l/Emr  » 


(4.14) 


and  thus  from  equations  (4.8,  (4*12)  and  (4.13) 


T  =  "  (Ek)m[\  (L)j^5«6  E,bl  (01  +  62)/Emr 


J 


-5.6  (0t  +  02 )bl/r 
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-  It-t  ({  T  z  (T  -  0  )  a  A  , 
9T  7  z  N  z  m7 


5.6  (9,  +  e.)  [(Ek)m  61  (Em  -  E,em)/Emr  +  ff£ 


'm  '  m 


j  TZ(L  -  L)  zcLaJ 


The  integral  expression  may  be  written 


f 


Tz  (L  -  L)  zcLA  =  -  6  Tzz2dA  , 


and  to  a  sufficient  approximation  the  term  Tg  in  this  way  may  be  replaced  by 

6  and  thus 
in 


T'=  5.6  (6i  +  02  )  ^  [jEk)m  E0/Em  + 


9Ek  ~| 

9T  6mJ 


¥e  have 


V  =  j 


z2Ez  dA  , 


and  to  a  sufficient  approximation  this  may  be  written 


d  =  Em!  , 


and  thus 


1/r  =  -  £  =  5-6  (e,  ♦  e2)  [(Ek)m  Eo/E^  +  ^  eB 


/eb  ,  (4-.15) 


dk 


and  this  may  be  written  (since  km  =  k^,  +  0m) 


V*  =  5.6  (e,  ♦  %)  |  [eb  *  |?  §  emJ 


(4.16) 


It  remains  to  calculate  the  position  of  the  reference  axis  z  =  o. 

To  simplify  the  solution  we  suppose  the  internal  contour  of  the  cylinder  is 
a  perfect  circle  of  radius  r.  Thus 

|zdA  =  j2  (r  sin  0  -  a)  rd0  j^L  +  (z  +  a)^j  , 


where  a  is  the  distance  of  the  reference  axis  z  =  o  above  the  diameter  of  the 
circle .  Thus 


/ 


z  d  A  =  b (irr2  +  2va2  )  -  2rrra  (L  +  —  a)  . 


Equating  this  to  the  value  given  by  equation  (4*14)  we  find  (neglecting  a2 
compared  with  r2 ) 


a 


br  (1  -  KL) 

t  r^icb2 


(4.17) 


where  K  =  5*6  (©i  +  02)E1/Em  . 

The  temperatures  in  2-inch  and  3“inch  motors  have  been  determined 
theoretically  (Ref.  1 )  and  the  temperatures  on  the  outside  #f  the  tube  have 
been  determined  experimentally  at  A.R.D.E.  and  P.D.E.  (Ref.  2).  The  measure¬ 
ments  and  calculations  have  been  confined  to  three  positions  along  the  tube 
namely  (i)  near  the  head  end  (ii)  near  the  nozzle  end  and  (iii)  at  the  central 
section.  At  the  nozzle  end  and  at  the  middle  sections  the  theory  and 
experiment  agree  reasonably  well  but  at  the  head  end  the  observed  outside 
temperatures  are  appreciably  greater  than  those  calculated.  The  theoretical 
values  of  the  temperature  at  the  inside  wall  are  readjusted  using  a  heat 
transmission  coefficient  which  gives  the  observed  outside  temperature.  With 
sufficient  accuracy  0  may  be  replaced  by  (01  +  02)/2.  The  radii  of  curva¬ 
ture  have  been  calculated  for  the  3“ inch  rocket  with  the  tubular  charge  at 
20  G  and  with  b  =  10”3  ins.  The  curvatures  are  given  in  the  following  table, 
R  being  in  inches. 


Nozzle  End  c  Head  End  c  Middle  Section  _ 

n  a  /n  a  r\J  cs  .  q  a  /n  v  AnP  a  o  A  /r>  ^  An? 


t  (secs.) 

0^  +  02 

l/R  x  10^ 

©i  +  02 

l/R  x  10^ 

.  ®2 

l/R  x  105 

0.05 

276 

1.31 

112 

0.49 

194 

0.88 

0.15 

608 

3.46 

292 

1.40 

488 

2.48 

0.25 

858 

5.50 

406 

2.08 

658 

3.68 

0.35 

982 

6.66 

482 

2.56 

792 

4.76 

0.45 

1108 

8.20 

554 

3.01 

900 

5.70 

0.55 

1200 

9.76 

610 

3.35 

988 

6.56 

0.65 

1284 

12.20 

672 

3.75 

1068 

7.51 

0.75 

1350 

X 

700 

3.93 

1148 

8.23 

0.85 

1398 

ii 

732 

4.15 

1228 

10.14 

0.95 

1440 

X 

762 

4.34 

1280 

11.87 

3iE0/Em  is  not  known  accurately  at  these  temperatures  -  it  increases 
very  rapidly  giving  rise  to  very  high  bending. 

At  a  specified  time  the  curvature  varies  along  the  length.  We  assume 
it  can  be  expressed  in  the  form  of  a  quadratic,  or 

"  =  1/R  =  a-,  +  b^  x  +  c^  x2 


Writing  1/R1  ,  1/r2  and  l/R3  as  the  curvature  at  the  head  end,  the  middle 
section  and  the  nozzle  end  respectively,  we  can  express  a1 ,  b1  and  c1  in  terms 
of  these.  The  angle  between  the  nozzle  axis  and  the  rocket  axis  at  the  head 
is  thus 

l 

=  a,  l  +  b,  l2/2  +  c1  l3/3  =  (1/Ri  +  4/R2  +  1/R3 )  1/6  , 

x  =  o 


pa 

L3*. 
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I  being  the  length  of  the  motor*  This  angle  is  calculated  at  the  various 
times 


t  (secs.) 

0.05 

0.15 

0.25 

0.35 

0.45 

0.55 

O.65 

0.75 

0.85 

0.95 

% 


Angle  between  shell  axis  and  nozzle  axis  due  to  bending 

0.41  x  10"^  radians 

1.11  "  " 

1.70  "  " 

2.12  "  " 

2.53  "  " 

2.95  "  " 

.  3.45  " 

3.70  "  "  * 

4.56  "  "  * 


5.34  " 


"  St 


*These  are  taken  to  be  l/R2»  this  underestimates  the  angle. 


In  the  bending  experiments  referred  to  in  other  sections  the  head  end 
of  the  rocket  is  rigidly  clamped  and  the  lateral  movement  of  the  venturi  end 
measured.  This  lateral  movement  is 


[y]0l  =  a,  l2/2  +  l3/6  +  Cl  1 4/l 2  =  l2  [1/6R,  +  1/3R*]  ,  (3.18) 


The  calculated  movements  for  the  3“inch  rocket  with  tubular  charge 
and  b  =  10“3  ins.  are  given  in  the  following  table. 

t  (secs.)  Lateral  movement  of  nozzle  end  with  head  end  clamped  (ins.) 


0.05 

0.15 

0.25 

0.35 

0.45 

0.55 

O.65 

0.75 

0.85 

0.95 


0.007 

0.021 

0.032 

0.041 

0.049 

0.056 

0.063 

0.069 

0.082 

0.095 


Due  to  the  curvature  of  the  tube  the  centre  of  gravity  moves  laterally. 
If  y  is  the  lateral  movement  of  the  C.  of  G. ,  dt  the  lateral  movement  of  the 
end  x  =  l  and  <f>  the  angle  between  the  shell  axis  and  the  nozzle  axis,  the 
effective  tolerance  angle  is 

a  =  <f>  -  (di  -  y)/l  -  e)  , 

e  being  the  distance  of  the  C.  of  G.  from  x  =  o.  The  tenn  y  is  given  by  the 
equation 

(W1  +  wl  +  W2  l2 )  y  =  w  [a,  l3/6  +  b,  l4/24  +  c,  l5/60] 

+  W2  [(a,  l  +  b,  l2/2  +  c,  l3/3)l2  +  (a,l2/2  +  b,  l3/6  +  c^4/^)] 


23 


the  notation  being  as  in  section  3*  The  calculated  values  are  given  in  the 
following  table. 

t  Equivalent  nozzle  offset  (radians  x  10 


0.05 

0.28 

0.15 

0.71 

0.25 

1*09 

0.35 

1.32 

0.45 

1.61 

0.55 

1.88 

0.65 

2.17 

0.75 

2.46 

0.85 

3.00 

0.95 

3.91 

Approximately  the  effective  tolerance  angle  is  a  linear  function  of 
the  time.  Further  the  velocity  of  the  rocket  is  approximately  a  linear 
function  of  the  time.  We  have 

a  =  /?  +  yt  ,  V  =  V0  +  ft  ,  f  =  F/W  . 

V0  being  the  launch  velocity. 

«  *  /?  +  y  (V  -  V0  )/f  -  A  +  y1  V  . 

The  yaw  at  a  distance  s  is 


if  fs  Min 

is.  x-v 


.)  •  / 
— 7  sm  n  (s 


-  u)  du  , 


which  may  be  written 


1  f  &  F  ( l  **  e)  .  f  s 

6  =  SV  /  Ii  V”  &  sin  n  (s  “  u)  du  +nv  /  Fy^1  “  e)  sin  n  (s  -  u)du. 
J  So  Jsn  Ii 


The  yaw  due  to  a  constant  cross-wind  w1  is 


6  =  ■^L  (l  -  cos  n  (s  -  s0  )} 


Thus  if  we  assume  the  thrust  to  he  constant,  the  yaw  due  to  the 
variable  equivalent  nozzle  offset  a  is  equal  to  the  yaw  due  to  a  constant 
nozzle  offset  /?*  added  to  the  yaw  due  to  a  constant  cross  wind  of  strength  w1 
where 

_  F.  yi  ( l  -  e) _ F_  (l  -  e)  z  _  W(l  -  e)y/ 

Wl  "I,n?  “I^2  f  -  1in 
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W  being  the  total  mass  of  the  rocket.  Thus  if  X,  is  the  deviation  of  the 
rocket  in  flight  due  to  a  nozzle  offset  of  10~3  radians  and  X2  that  to  an 
unit  cross-wind,  the  deviation  due  to  the  variable  nozzle  offset  a  is 

5  =  (/?  -  yV0/f )  X,  x  103  +  W  (l  -  e)  yXa/^n2  . 

For  the  3-inch  rocket  with  parallel  sided  shell,  A.D.  fuze  and  tubular 
charge  the  values  of  X1  and  X2  are  0.8  degs.  and  0.12  degs.  respectively. 

Thus  the  deviation  of  this  rocket  with  wall  thickness  variation  specified  by 
b  =  10”*3  ins.  is  1.0  degs. 

In  all  the  above  work  it  has  been  assumed  that  the  elastic  limit  is  not 
exceeded,  that  is,  Hooke's  law  is  valid  apart  from  the  variation  of  Young's 
modulus  with  temperature .  The  problem  of  thermal  stresses  in  3-inch  and 
2-inch  rocket  tubes  has  been  dealt  with  in  detail  by  Oldroyd  (Ref.  3) •  From 
his  analysis  he  concludes  that  in  firings  at  air  temperature,  serious  excesses 
of  the  calculated  stresses  above  the  elastic  limit  occur  near  the  surfaces  of 
the  tube  wall  but  the  major  part  of  the  maeterial  of  the  tube  is  not  stressed 
beyond  this  limit.  The  highest  stresses  occur  towards  the  venturi  end  of  the 
combustion  chamber  and  are  compressive  stresses  near  the  inner  surface  and 
tensile  stresses  near  the  outer  surface  of  the  tube  wall.  Beyond  the  elastic 
limit  the  value  of  Young's  modulus  decreases  very  rapidly  with  increasing 
strain  and  also  we  may  expect  that  the  coefficient  of  thermal  expansion 
increases  rapidly.  The  expression  for  T  given  by  equation  3*4  is  a  measure 
of  the  variability  of  kET  around  the  section.  If  we  exceed  the  elastic  limit 
then  clearly  very  much  larger  variations  of  kE  are  possible  and  this  we  may 
expect  to  have  a  very  serious  effect  on  bending.  The  stress  system  is 
exceedingly  complex  and  no  theoretical  calculations  are  possible;  in  any  case 
the  variability  depends  to  a  marked  extent  on  the  residual  stresses  in  the  tube. 

The  effect  of  thermal  stresses  is  best  illustrated  by  comparing  tubes 
of  different  wall  thickness.  The  standard  3-inch  tubes  considered  by 
Oldroyd  are  0.080  ins.  thick.  We  consider  also  tubes  of  the  same  internal 
diameter  but  with  wall  thickness  0.20  ins.  The  temperatures  in  the  thick 
tubes  have  been  calculated  using  the  same  emissivity  coefficients  as  used  in 
calculating  the  temperatures  of  the  standard  tubes  so  that  the  results  are 
strictly  comparative.  The  temperatures  are  calculated  for  uncoated  tubes 
with  tubular  charges  and  the  venturi  end  of  the  combustion  chamber. 

Tube  Temperatures  (°C.) 

0.08  in.  tube  0.20  in.  tube 


t  (secs.) 

Inner  surface 

Outer  surface 

Inner  surface 

Outer  surface 

0.2 

520 

260 

490 

8 

0.4 

660 

460 

580 

70 

0.6 

760 

580 

630 

150 

0.8 

830 

670 

- 

- 

The thermal  stresses  are  calculated  using  a  simplified  version  of 
Oldroyd' s  treatment  for  both  wall  thicknesses.  We  make  two  assumptions,  (i) 
there  is  cylindrical  symmetry  about  the  rocket  axis,  (ii)  if  v  is  the  displace¬ 
ment  in  the  axial  direction,  we  assume  is  constant  through  the  wall  of  the 
tube.  Both  these  conditions  are  approximately  time. 

If  u  is  the  displacement  in  the  radial  direction  i.e.  an  element 
originally  at  a  radius  r  is  displaced  to  r  +  u,  the  classical  equations  of 
elasticity  take  the  form 


25 


> 


(  _rrr 


*■#2  [g*f  ♦<--'>£-<.♦  n»«]  • 


G  being  the  shearing  modulus,  m  the  inverse  Poisson's  ratio  and  R,  H  and  X 
are  the  radial,  hoop  and  axial  stresses  respectively.  Prom  the  condition 
of  equilibrium  between  these  stresses  we  have 


H-R  =  rfS  . 


Prom  this  condition  we  find 


JL  (1  JL  (ur))  -  5LLl 

3r  vr  3r  v  n  m  -  1  3r 


On  integration  we  have 


—  -  (ur) 

r  3r  N  ' 


ur 


=  2-Ul  f, 

m  -  1  I 


r  k  T  d  r  +  C0r*/2  +  Ct  . 


The  constants  C0  and  C1  are  obtained  from  the  conditions  that  R  =  -  P 
(the  internal  pressure)  at  the  inner  radius  and  R  =  0  at  the  outer  radius  and 
thus  we  obtain  the  expressions  for  R,  H  and  X. 

The  thermal  stresses  given  in  the  following  tables  are  calculated  on 
the  assumption  that  Hooke's  law  is  valid,  it  being  assumed  that  the  elastic 
limit  is  not  exceeded,  the  values  of  Young 1 s  modulus  being  appropriate  to  the 
temperature  of  each  element. 

Hoop  stresses  at  nozzle  end  of  the  tube  (tons  per  sq.in) 


t 


(secs. ) 

0.08  in.  tube 

Inner  surface  Outer  surface 

0.20 

Inner  surface 

in.  tube 

Outer  surface 

0.2 

-52.0 

51.0 

-100 

32.0 

0.4 

-40.4 

26.1 

-108 

45.1 

0.6 

-28.0 

20.3 

-107 

48.7 

+  sign  denotes  Tension,  —  sign  denotes  Compression. 

The  following  table  gives  the  distribution  of  stress  through  the 
wall  of  the  0.20  in.  tube  at  t  =  0.4  secs. 


26, 


Distance  from 
inner  surface  (ins. 

Hoop  Stress 
)  (tons/sq.in. ) 

Longitudinal  Stress 
(tons/sq.in. ) 

0 

-108 

-110.0 

0.02 

-70.9 

-73-5 

0.04 

-45.1 

-42.7 

0.06 

-15.4 

-18.0 

0.08 

+4.1 

+  1.5 

0.10 

+18.9 

+16.3 

0.12 

+53.4 

+30.8 

0.14 

+37.2 

+34.6 

0.16 

+41.9 

+39.3 

0.18 

+44.4 

+41.8 

0.20 

+45.1 

+42.6 

The  following  table 
normal  temperature  s . 

gives  the  strength  properties  of  the  rocket 

Longitudinal  (t.s.i.) 

Transverse  (t.s.i 

Elastic  limit 

26 

21 

0.1$  Proof  Stress 

38 

37 

0.2$  Proof  Stress 

40 

40 

U.T.S. 

44 

44 

It  is  seen  that  the  calculated  stresses  are  markedly  in  excess  of  the 
elastic  limit  with  "both  0.08  in.  and  0.20  in.  tubes  and  clearly  with  both 
tubes  there  is  plastic  straining,  the  effect  with  the  0.20  in.  tube  being  very 
markedly  greater.  The  figures  apply  at  the  venturi  end  of  the  tube.  At  the 
middle  section  of  the  tube  the  calculated  stresses  are  approximately  half 
those  given  above  and  with  0.08  in.  tube  there  is  only  slight  plastic  deforma¬ 
tion  at  this  section.  With  the  0.20  in.  tube  the  plastic  deformation  is  still 
considerable  at  the  middle  section.  With  coated  0.08  in.  tubes  the  calculated 
stresses  exceed  the  elastic  limit  so  an  appreciable  extent  only  near  the 
venturi  end  but  with  the  coated  0.2  in.  tube  the  elastic  limit  is  exceeded 
appreciably  beyond  the  middle  section  of  the  tube. 


APPENDIX  I 


Calculation  of  the  heating  of  the  rocket  tube  in  the  case  of  star  centre 

charges  which  are  sealed  at  the  head  end 

We  make  the  assumption  that  on  ignition,  propellant  gases  flow  round 
the  charge  to  pressurise  the  annular  space  between  charge  and  tube.  Plow  is 
assumed  to  cease  and  then  interchange  of  heat  between  the  stagnant  gases  and 
the  tube  wall  tubes  place.  The  charge  coating  is  assumed  to  absorb  negligible 
heat. 


The  following  constants  are  taken  for  the  gas 


Pressure  =  70  atmospheres 

Universal  gas  constant,  R,  =  0.08204  litre  atmospheres/  C. 

Gram  molecules/grm.  n  =  0.0424  moles/grm. 

Specific  heat  of  gas  at  constant  pressure,  C  =  0.369  cals/grm. 


Isobaric  flame  temperature  degrees  above  ambient  T^  =  2400 


centigrade  degrees 

Also;-  The  gap  between  propellant  coating  and  wall  d  =0.2  cms. 
Density  of  steel,  q  =  7*85  grm/cc. 

Specific  heat  of  steel,  s  =  0.12  cals/grm. 

Thickness  of  steel  tube,  t  =  0.2  cms. 


We  calculate  the  specific  volume,  V, 


nRTp  .0424  x  .08204  x  2400  litre/grm. 
P  =  70 


V  =  120  cc/grm. 

d.  c  T 

Heat  content  of  the  gas  per  square  cm.  of  surface  =  p  p 

V 

0.2  x  0.369  x  2400  „  ,  -,r  .  , 

- - -  =  1.476  calories/sq.cm 

If  the  equilibrium  temperature  is  T,  then  we  have 


dCT  dCT  .  _ 

p  p  =  p  +  t  q  s  T 


V 


V 


°'-2,,X  °‘?20  X  2400  -  °-240-^  +  °-2  *  7-®5  *  °-12  T 

from  which  T  =  7»8°C. 

In  fact,  the  flow  will  cause  the  tube  to  heat  to  a  higher  temperature 
at  the  source  end,  and  lower  at  the  sealed  end,  but  the  average  temperature 
rise  along  the  tube  should  be  about  7.8  C.  Some  complications  occult  due  to 
gas  dissolving  in  the  charge  restrictive  coating  causing  more  gas  to  heat  the 
tube  than  calculated  above.  The  distance  d  varies  round  the  circumference 
and  some  asymmetric  heating  arises. 
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APPENDIX  II 


The  effect  of  variations  of  the  wall  thickness  and  coating  thickness  on  the 

temperature  of  the  tube 

It  can  be  shown  that  the  transfer  of  heat  through  the  tube  during  the 
burning  period  is  almost  entirely  radial  and  there  is  negligible  conduction 
around  the  circumference.  The  wall  thickness  is  small  compared  with  the 
radius  of  the  tube  and  thus  the  problem  is  virtually  one  dimensional. 


steel 


X»«2 

We  use  the  following  notation 
0  is  the  local  temperature 

t  =  o  is  the  interface  between  the  coating  at  the  gas, 
i  =  is  the  interface  between  the  coating  and  the  steel  tube, 
i  =  &2  is  the  outer  surface  of  the  tube. 


®i  »  Pi  >  k,  ,  specific  heat,  density  and  thermal  conductivity  of  the  coating 
s2,  p2,  k2 ,  specific  heat,  density  and  thermal  conductivity  of  the  steel. 


K, 


,  Kj  = 

Si  Pi  S2 p2 


Z  is  the  temperature  of  the  gas. 

The  equation  of  heat  conduction  in  the  coating  is 


0i 


01  1  0t 


and  in  the  steel 


ax  ^  ax  )  =  (S2P20) 


0X 


0X 


0t 


During  the  period  of  burning  there  is  negligible  heat  transfer  from  the  outer 
surface  of  the  steel  thus 


96  „ 

0T  =  0  »  1  = 


(i) 


thus 


At  x  =  a.A  the  temperatiire  and  of  course  the  heat  flow  are  continuous 


6  =0 
a5  -  e  a1  +  e 


(ii) 
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.  /96\ 

-  (gi) 


a<  +  e 


(iii) 


e  "being  a  very  small  positive  quantity. 

At  1=0  there  is  heat  transfer  from  the  hot  gases  to  the  coating  and 
this  gives  the  boundary  condition 

t,  g  =  m(e-»)  .«  =0  , 

where  h^  is  the  transmission  coefficient  from  the  hot  gases*  This  is  usually 
written 

=  h  (0  -  z)  ,  x  =  o  ,  (iv) 

dx 

The  quantities  k,  s  and  p  depend  on  the  local  temperature  and  the 
value  of  depends  on  the  boundary  temperature  and  the  looal  gas  velocity  and 
these  quantities  vary  appreciably  during  burning.  The  temperatures  in 
standard  tubes  have  been  obtained  accurately  by  numerical  methods.  The 
amount  of  work  involved  is  immense  and  since  we  are  concerned  with  the  effects 
due  to  small  changes  in  parameters,  appreciable  errors  may  arise  if  we  use 
these  methods.  Thus  an  alternative  method  has  been  adopted. 

In  the  first  instance  we  assign  mean  values  to  s,  p,  k  and  h.  The 
differential  equations  thus  reduce  to 


a20  _  i  36 
te*  =  K,  3t  * 


020  1  30 

9I7-  =  K,  0 1  ’ 


o  <  i  <  a., 


a1  <  x  < 


(v) 


(vi) 


The  solution  of  equations  (v)  and  (vi)  subject  to  the  boundary 
conditions  (i)  and  iv)  is 


0  =  z  (1  -  A,  e  a'  -  Aj  ^  e  a2t 


(vii) 


where 


<f>  =  sin  Knx  +  an  cos  anx 

n  7!c7  h75q  7ST  '  <  1  <  1  ’ 


(sin  a  a,  a  cos  a  a,  )  cos  **  =0 

7§7  +  hVl,  757  TkT 


cos  an(a2  -  ) 


a,  <x<&z 


“TkT 


and  the  terms  a  are  the  roots  of  the  determinant 
n 
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=  0 


sin  aa,  a  cos  a  a,  cos  a  (a*  -  a, ) 

7k,  iTk,  VK,  *  7^ 


k,  ,-cos  aa,  _  a  sin  aa  k^  sin  _a_  (a2  -  a,  ) 


■Sk.\  ^  7k,'  hTK,  7k,  * 


k 


7h 


Also 


and 


^^TVk,  A3'5’  f  '  *n  +  Sll>‘  f  V  dl) 

n  /  J0  Jai 

p  4«4c  .[>  3in  ^n8’  +  an  (*  +  &  3in  ^n1’) 

I  n  L2  tan  TET  D?  2  ’“n  "7E7 


i  / .  uua  2a1a1  \  | 
+  2h(l~  ^  >  J 


f  **  .  2  ,  (sin  a  a,  ag  cos  aa,  \2 
ia,  *”  "  =  TC^Pt, 

|~cos  g(a^  -  a, ) 


fa,  -  a,  VKg  sin  2g(a;>  -  a,  )~| 

L  2  +1*n  J 


In  all  these  expressions  an  is  denoted  by  g 


but  with 


In  the  absence  of  the  coating  the  solution  is  given  by  equation  (vii) 


6  =  sin  ax  +  a 

rn  n  __n 

h 


cos  a2x 


+  an  )(aa2 
h5. 


+  sin  ana, 


cos  a2a) 


In  the  first  instance  we  take 

S,  S  0*2  P,  2  1*9^  lc,  2  0*0036 

S2  =  0.135  Pz  =7.7  a  0.098 

these  being  in  C.G.S  and  °C  units.  The  mean  thickness  of  the  coating  is 
0.006  ins.  and  the  mean  wall  thickness  0.08  ins. 

Three  values  of  h  for  the  coated  tubes  have  been  used  in  the  calcula¬ 
tion  namely  h  2  10,  25  and  75  and  these  cover  the  range  of  h  in  practise. 
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7e  proceed  thus.  For  each  value  of  h  we  calculate  the  mean  tube 

temperature  for  three  values  of  aj  -  a,  ,  namely  the  value  0.203  cms  (0.08  ins.) 

and  two  other  values  0.183  cms  and  0.220  cms.  Let  (0  )  and  (0  )  he  these 

m  1  m  2 

two  latter  values  and  L  the  wall  thickness,  thus 


(90  2\ 
'■0L  ; 


L  =  0.203 


(em>2  -  (8m>, 


J  / (0.220  -  0.185) 


In  addition  we  calculate  01  and  02  the  temperatures  at  x  =  a1  and 
i  -  a2  with  ag  -  a1  =  0.203.  It  is  found  that  at  the  three  values  of  h  the 
term 


is  practically  independent  of  the  time  and  of  h. 


Coated  tubes 

X 


t  (secs.) 

h  =  10 

h  =  25 

h  =  75 

0.15 

3.97  “ 

3.97  * 

3.97  ! 

0.25 

4*34 

4.35 

4.35 

0.35 

4.50 

4.47 

4.45 

0.45 

4.54 

4.48 

4.36 

0.55 

4.55 

4.46 

4*34 

0.65 

4*54 

4.42 

4.22 

0.75 

4.55 

4.35 

4.14 

0.85 

4.50 

4.34 

4.04 

0.95 

4.45 

4.26 

3.94 

1.05 

4.39 

4.21 

3.83 

"The  series  expressions  converge  slowly  at  t  =  0.15  and  these  may  be  in  slight 
error. 

The  calculations  have  been  repeated  with  the  uncoated  tube  with  h  =  0.5, 
h  =  1 ,  h  =  2,  these  covering  the  ranges  of  h  obtained  in  practise. 

Uncoated  tubes 
X 


(secs. ) 

h  =  0.5 

h  =  1.0 

CM 

II 

0.15 

3.95 

4.12 

4.16 

0.25 

4.23 

4.30 

4.34 

0.35 

4.44 

4.36 

4*34 

0.45 

4.49 

4.38 

4.24 

0.55 

4.47 

4.33 

4.13 

0.75 

4.42 

4.22 

3.84 

0.95 

4.38 

4.08 

3.55 

1.05 

4*34 

3.92 

3.26 
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We  are  concerned  with  the  first  0*6  secs,  of  burning  and  in  this 
region  X  has  a  constant  value  for  both  coated  and  uncoated  tubes.  The  reason 
for  this  constancy  of  X  is  apparent.  Let  us  suppose  that  amount  of  heat 

crossing  the  boundary  into  the  tube  is  independent  of  small  variations  of  the 
wall  thickness.  Thus  the  mean  temperature  of  the  tube  is 

e  H 

m  Ls2  pz 

and  if  we  assume  s2p2  as  constant  where 


9L  l2s2Pz  “  I*  m 


or 


1 

fl  — m  = 
m  3L 


1_ 

L 


We  have  L  =  0.203  so  that  1/L  =  4*9»  The  value  (01  +  02)/2  is  approximately 
equal  to  0m  and  thus  we  see  that  the  constancy  of  X  is  due  to  the  condition 
that  the  heat  that  has  previously  crossed  the  boundary  varies  only  slightly 
with  small  variations  in  the  wall  thickness.  The  temperature  01  and  02  have 
been  used  as  reference  temperatures  since  these  have  already  been  computed 
accurately. 

Repeat  calculations  have  shown  that  X  is  invariant  to  small  changes 
in  k,  and  as  we  may  expect  and  thus  we  can  confidently  use  the  formula 

®®m  _  X( 0<  +  0g ) 

0L  2 

and  the  value  of  X  taken  is  4.4  in  centimetre  units  or  11.2  when  the  wall 
thickness  is  measured  in  inches.  This  is  correct  for  both  coated  and 
uncoated  tubes. 

The  values  of  0.)  and  02  have  been  computed  accurately  in  earlier  works 
(Ref.  1 )  and  the  results  are  given  in  Table  I. 
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APPENDIX  III 


Bending  of  the  motor  in  static  firing 

In  Section  2  we  consider  the  "bending  of  the  motor  tube  in  flight. 

There  it  is  noted  that  the  natural  frequency  of  bending  is  87  cycles  per 
second.  In  the  static  firing  the  head  end  of  the  motor  tube  is  rigidly 
clamped  and  this  changes  the  natural  frequency.  We  have  seen  in  Section  2 
that  for  a  perfect  rocket,  the  pressure  enters  only  in  the  gas  whip  terms  and 
provided  we  take  the  values  of  Young's  modulus  appropriate  to  normal  tempera¬ 
tures,  this  affects  the  amplitude  terms  having  only  a  secondary  effect  on  the 
frequency.  Further  the  thrust  has  only  a  small  effect  on  the  frequency  under 
these  conditions.  The  detailed  problem  is  considered  in  a  later  report. 
Neglecting  the  pressure  and  thrust  terms  the  differential  equation  becomes 


Dy 


+  X2 


¥- 


W 


X  =  EI 


D=£ 


(D 


w  being  the  weight  per  unit  length  of  the  shaft.  With  sufficient  accuracy 
we  may  neglect  the  length  of  the  venturi  and  assume  the  venturi  assembly  to 
be  a  point  mass  acting  at  its  centre  of  gravity.  The  boundary  condition  thus 
become 


EI  D3y  =  W  |^f  ,  I  =  l 


(2) 


Dy  =  o 


,  X  =  i 


(3) 


W  being  the  mass  of  the  venturi  assembly  and  l  the  distance  of  its  centre  of 
gravity  from  the  clamp.  In  addition  the  boundary  conditions  at  the  clamp  are 


y  =  Dy  =  o 


x  =  o 


ipt 


Writing  y  =  Ye  ,  the  general  solution  of  the  differential  equation 
satisfying  the  conditions  (4)  is 


Y  =  A 


Sinh  z1  -  sin  z1 

+  B  Cosh  z1  -  cos  z1 

_  — . 

L 

where 


zi  =  (pX)^x 


Substituting  this  into  the  boundary  conditions  (2)  and  (3)  and  using  the 
notation 

S  =  Sinh  z  ,  s  =  sin  z  ,  C  =  Cosh  z  ,  c  =  cos  z  ,  z  =  (pX)^"l 

we  obtain  the  equations. 

The  determinant  of  these  equations  reduces  to 


1  +  cC  +  (Sc  -  Cs)  =  o 


(5) 


W 


This  is  an  equation  in  z  and  W/wl.  For  a  given  value  of  z  we  calculate  — 
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The  practically  important  values  are  given  in  the  following  table 


wl 

z 

W 

0.2 

0.00053 

0.4 

0.00854 

0.6 

0.0437 

0.8 

0.1411 

1.0 

0.3620 

1.2 

0.8263 

1.3 

1.230 

1.4 

1.841 

1.5 

2.833 

1.6 

4.57 

1.8 

22.62 

For  a  specified  value  of  —  we  determine  z  either  by  interpolation  or 
graphically.  The  natural  frequency  is  thus 

1  z2  ,ElJ 

a?  F  * 


1 
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APPENDIX  IV 


The  stability  of  the  rocket  charge 

The  rocket  charge  rests  on  the  grid  and  usually  it  is  not  anchored  at 
the  head  end  "but  depends  on  the  tube  for  support.  In  the  case  of  tubular 
charges,  plastic  tabs  placed  on  the  charges  restrict  the  independent  motion. 

It  is  of  interest  to  determine  the  stability  of  the  charge  when  it  is  supported 
only  at  one  end.  Denoting  the  Young's  modulus  of  the  charge  by  E  and  the 
moment  of  area  by  I,  the  bending  moment  is 


The  longitudinal  force  is 


Gjc  =  PqAo  +  wxf 

P0  being  the  stagnation  pressure,  Aq  the  cross-sectional  area  and  w  the  mass 
per  unit  length  of  the  charge.  The  shearing  stress  is  given  by 

0  3B  „  3y 
s  +  3^  +  g>  =  0  ’ 

and  in  the  absence  of  external  forces  and  couples  the  equation  of  motion  is 


9  y  wfx  92y  wf  9y  w  92y 

9x2  li  9^  +  ti  9  x  + 


=  o 


Writing  •  —=  a.  ,  ^  =  c2  we  get 


D4y  +  aD(xDy)  +  f^r  =  o 

Let  us  suppose  y  contains  t  in  the  form  e^c^;  we  have 


D4y  +  aD(xDy)  +  \2y  =  o 

In  the  first  instance  let  us  suppose  the  nozzle  end  of  the  charge  to 
be  rigidly  connected  to  the  tube,  while  the  head  end  is  free.  Thus 

y  =  U  =  o  ,  x  =  l  ;  B  =  o,S  +  G-0|^  =  oatx=o; 

These  may  be  written 

,.§..*..1,  =  o  at  *  =  l. 

The  solution  for  y  may  be  expressed  in  the  series  form 


y  =  S  An,n/n! 
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Substituting  this  into  the  differential  equation  gives  the  recurrence 
relations  for  the  terms  A  .  All  the  coefficients  may  be  expressed  in  terns 
of  the  first  four.  The  'Boundary  conditions  then  give  four  relations  which 
Ao  »  A(  ,  A2  and  A3  must  satisfy.  The  determinant  of  these  equations  may  be 
expanded  in  the  form 


+  X4 


+  X6 


al3 

“irH 

4 

L Zl 

(al3)3  -ff 

(al3) 

3  ♦  f§  (al3 

'Y  (al3) 

5  . 

(*- 

18  al3  .  200  (al3)2_. 

7,  at  +  10J  v.at  ) 

(a.3/ . 

[1- 

192 
11 1 

ai3  +  *** 

at  +  141 

(al3 )2 

2fl_ 

4. 

JJ2J 

1  •  1 

+ 


If  the  numerical  values  of  a  and  l  are  substituted,  the  equation  A  =  0 
gives  us  the  values  of  X2.  If  all  these  are  negative,  the  charge  is  stable. 

If  however  one  or  more  positive  values  of  X2  are  obtained,  the  charge  is  unstable. 
The  first  root  X2  may  be  obtained  quickly  from  the  following  approximate 
equation 

X?  I4  =  -  12.37  +  1.58  al3  . 


For  the  3-inch  rocket  No.  1  tubular  charge  we  have 

f  =  2000  f.s.2 
w  =  3.54  lb/ft., 

Q 

E  =  1.6  x  10  pdls. /sq.ft. 

I  =  1.25  x  104  ft.4  , 

l  =  3.6  ft. , 

thus  al3=  16.5 

The  greatest  stable  length  is  such  that  X<  =  o  so  that  it  is  given  by 

1.58  a  l3  =  12.37  or  a  l3  =  7.84 


thus  the  charge  is  completely  unstable .  The  length  of  the  charge  that  is 
just  stable  is  2.8  ft. 

Similarly  if  we  supposed  that  in  addition  the  head  end  of  the  charge 
were  pin  pointed  so  that  y  =  B  =  o  at  x  =o,we  find  the  charge  is  still 
unstable. 


If  the  charge  were  rigidly  located  at  the  head  end  so  that  y  = 

O* 

a  0  at  1  =  o  and  x  =  l ,  the  charge  is  almost  neutrally  stable . 
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FIG.  3 


*S1IW  '13SJ40  31ZZON  1N31VAM03 


I— 

o 

o 


S3H0NI  ‘ON 9  31ZZ0N  JO  JLN3W3A0W  1VH31V1 
CONFIDENTIAL 


FIG.  3.  BENDING  OF  3  INCH  ROCKET  ASSUMING  EVEN  HEATfNG  FROM  TUBULAR  CHARGE  &  WALL  THICKNESS 
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FIG.  4.  3  INCH  N  o.  I.  2-75-  0-75  SUK  DRILLED  TUBULAR,  -008  INCH  REFRACTORY  LINING  INTACT 
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FIG.  7 
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FIG.  7.  3  INCH  N«.l.  XII  SUK,  008  INCH  REFRACTORY  LINING  INTACT 
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FIG.  8.  3  INCH  ROCKET,  COATED  TUBE,  FILLED  SU/K/X/II. 
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FIG. 9 
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FIG.  9  3  ROCKET.  NO  INSULATION.  CHARGE  XII  SUK. 

TEMPERATURES  ON  OUTSIDE  OF  TUBE. 
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TEMPERATURES  ON  OUTSIDE  OF  TUBE. 
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FIG.  II. 
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FIG. 12.  3  INCH  No. I  XII  WITH  50  THOU  REMOVED  FROM  EACH  ARM  (jN  SUK)  50  THOU  BI02  NEOPRENE  LINING. 
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FIG. 14 
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FIG. 14.  MAXIMUM  TEMPERATURE  DISTRIBUTION  ON  3*'n?I  ROCKETS 
WITH  VARIOUS  LININGS  AND  XII  TYPE  CHARGES. 
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FIG.  15 
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FIG. 15.  3*N?I  ROCKET.  Xll-SUK  50  THOU  PLANED  OF  EACH  ARM,  BYI6B  SEAMLESS  LINING  25-60  THOU. 
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FIG. 16 
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FIG  I6.  MAXIMUM  TEMPERATURES  REACHED  ON  100  MM.  S.S.  ROCKET  FILLED  PROPELLANT 

F.565/I4/K  STAR  CENTRE  ^CHARGE  DURING  STATIC  FIRING. 

(AS  INDICATED  8Y  “tEM  PI  LAO*  PAINTS^) 
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FIG.I7 
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FIG.  17.  5"  No.l.  ROCKET  WITH  4315**  2-9*  F.B.2  /  S  U  K  STAR  WITH  OPEN  HEAD  END.  DISTRIBUTION  OF  MAXIMUM 
TEMPERATURES  REACHED  DURING  STATIC  FIRING  AS  INDICATED  BY  *TEMPILACT  PAINTS. 
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FIG.  18 
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FIG. 18.  3*R0CK£T  NO  INSULATION. CHARGE  FB2  STAR  SDK 
PLAIN  OBTURATOR.  HEAD  UNSEALED. 
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FIG. 19. 
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FIG.I9.  3*R0CKET  NO  INSULATION, CHARGE  FB2  STAR(SUK) 
PLAIN  OBTURATOR  HEAD  UNSEALED. 
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FIG.20 


6 

i 


CONFIDENTIAL 


FIG.20,  3  INCH  N9|  2-9  FB2  STAB  IN  SUK  REFRACTORY  LIN IN6  INTACT  AND  HEAD  END  OF  CHARGE  SEALED. 

fPLAIN  OBTURATOR.) 
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FIG.  21. 
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FIG.22. 
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BENDING  IN  INCHES. 
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END  OF  TUBE.  REFRACTORY  LINING  INTACT. 
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FIG.  24 
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FIG.  2  5 
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FIG. 2 5.  3 *BOCKET  CHARGE  F488/649  FB2  STAR  WITH  STEPPED  CONDUIT.  HEAD  END  OF  CHARGE  SEALED. 

HEAD  END  OF  TUBE  SEALED  BY  WELD. 
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FIG.26. 
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HEAD  END  OBTURATOR  WITH  CASHEW  NUT  POLYMER  /PARAFORMALDEHYDE  (RD  CEMENT  1286.) 
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FIG.  27. 
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FIG.  27  3  INCH  Nt.  I  2-75-  0  75  SUK  DRILLED  TUBULAR.  -008  REFRACTORY  LINING  INTACT.  *18  INCH  THICK  TUBES. 
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FIG.  28 
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FIG.  28.  3  INCH  No.l. 2-75- 0-75  SUK.  DRILLED  TUBULAR.  008  REFRACTORY  LINING  REMOVED.  •  18  INCH  THICK  TUBES. 
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FIG.  2  9. 
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•00*  REFRACTORY  LINING  REMOVED.  0-18  INCH  WALL  TUBE. 
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FIG.30. 
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FIG. 31. 


FIG.  31.  MEASUREMENT  Of  BENDING  OF  ROCKET  DURING  STATIC 

FIRING  BY  MEANS  OF  ‘FLOATING"  CAPACITY  GAUGE. 
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FIG.  32 
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FIG.  32.  MODIFIED  FIELDEN  UNIT  WITH  FLOATING  CAPACITY  GAUGE 
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FIG. 33 


FIG. 33.  TYPICAL  RECORDS  USING  MAGSLIP  TECHNIQUE. 

UPPER  TRACE  50  C.P.S.  TIMING. 

LOWER  TRACE  0  25  INCH  CALIBRATION  STEPS 
(FIRING  RECORD). 
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FIG.  34. 


FIG.  34.  MEASUREMENT  OF  BENDING  OF  ROCKET  OURING  STATIC  FIRING 

BY  MEANS  OF  EARTHED  CAPACITY  GAUGE. 
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FIG.  35. 
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ROCKET  MOVEMENT 
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FIG.  35.  MEASUREMENT  OF  BENDING  OF  ROCKET  DURING  STATIC 

FIRING  BY  MEANS  OF  MAGSLIPS. 
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